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INTRODUCTION 


The  objectives  of  the  SRI  hurricane  tracking  and  monitoring  re- 
search program  for  the  U.S.  Air  Force  Office  of  Scientific  Research 
were  to  collect  and  analyze  at  the  Wide  Aperture  Research  Facility  (WARF) 
the  HF  skywave  sea  backscatter  for  Gulf  of  Mexico  tropical  storms  and 
hurricanes,  and  to  determine  the  feasibility  of: 

• Tracking  tropical  storms  and  hurricanes  from  a long  distance. 

• Estimating  the  ocean  surface-wind  direction  field- 

• Estimating  the  significant  wave  height. 

• Estimating  the  ocean  surface-wind  speed. 

This  final  report  briefly  summarizes  the  research  program  and  results 
from  14  July  1976  to  30  November  1978. 

The  principal  results  of  this  research  are  summarized  in  more  de- 
tail in  the  technical  papers  submitted  or  accepted  for  publication. 

These  papers  are  given  in  Appendices  A through  E.  In  Appendix  A we 
describe  the  results  of  the  Eloise  experiments  in  which  we  demonstrated 
that  we  could  map  the  surface  wind  direction  field  around  a hurricane 
and  determine  the  hurricane  position  from  this  map.  In  Appendices  B 
and  C we  describe  the  results  of  the  tracking  experiments  conducted  for 
hurricane  Anita  and  tropical  storm  Debra.  In  Appendix  D we  describe 
the  method  and  results  of  estimating  wind  speed  and  significant  wave 
height  for  hurricane  Anita  and  tropical  storm  Babe.  In  Appendix  E we 
describe  the  method  and  results  of  measuring  significant  wave  height 
and  the  wave  frequency  spectrum  for  low,  moderate  and  high  wind  speed 
cases.  We  have  or  will  be  presenting  the  results  of  this  research  at 
technical  sessions  of  several  professional  society  meetings.*  ^ 

g 

Technical  Report  Number  One  is  a summary  of  preliminary  interpre- 
tation of  the  WARF  hurricane  measurements,  theoretical  simulations  of 
the  sea-echo  Doppler  spectra,  methods  of  computing  wind  velocity  and 
wave  height  from  the  sea-echo,  theoretical  simulations  of  the  effect  of 
scattering  patch  size  on  the  measurements,  and  the  effects  of  the  iono- 
sphere on  the  interpretation  of  the  data. 


In  September  1975,  we  mapped  the  surface  wind  direction  field  and 
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estimated  the  storm  center  of  Gulf  of  Mexico  hurricane  Eloise  using 

the  WARF  high-frequency  (HF)  skywave  radar. ^ We  compared  the  radar 
wind  directions  to  available  in  situ  wind  directions  measured  at  the 
National  Oceanic  and  Atmospheric  Administration  (NOAA)  National  Data 
Buoy  Office  (NDBO)  data  buoys  EB10  and  EB04^  and  found  agreement  to 
within  10°.  We  compared  the  WARF  radar  position  estimate  to  the  official 
National  Hurricane  Center  (NHC)  track  and  found  agreement  to  within 
35  km.  These  early  experiments  suggested  that  the  WARF  HF  skywave  radar 
might  provide  a new,  independent  source  of  surface  hurricane  data  that 
were  not  routinely  available  through  the  use  of  more  conventional  sen- 
sors. In  addition,  the  radar  was  capable  of  continuously  recording 
data  at  locations  specified  in  real-time  by  the  radar  operator.  Our 
subsequent  work  has  demonstrated  the  feasibility  of  estimating  the 
ocean  wave  frequency  spectrum,  ocean  surface  current  velocity,  and  sur- 
face wind  velocity  from  the  WARF  measured  sea  backscatter  from  ocean 
gravity  waves.  In  Sections  II-V  we  summarize  the  principal  results  of 
research,  the  WARF  HF  skywave  radar,  the  sea-echo  Doppler  spectrum, 
and  the  radar  measurements,  respectively. 


II  WARF  MEASUREMENT  CAPABILITY 


During  the  past  26  months,  we  collected  and  analyzed  data  for 
hurricanes  Anita  (1977),  Babe  (1977),  Kate  (1976),  tropical  storm  Debra 

(1978),  and  further  analyzed  the  data  for  Eloise  (1975).  We  developed 

12  13 

tracks  for  Anita  (Appendix  B),  an  intense  hurricane,  and  for  Debra 

(Appendix  C) , a weak  tropical  storm,  from  the  WARF-derived  surface- 
wind  direction  maps.  The  mean  difference  between  the  WARF-determined 
positions  and  the  National  Hurricane  Center  (NHC)  interpolated  posi- 
tions along  the  NHC  track  was  less  than  20  km. 

12 

The  WARF  surface-wind  direction  maps  made  for  hurricanes  Anita 
9 

(Appendix  B)  and  Eloise  (Appendix  A)  were  compared  to  wind  direction 

maps  computed  from  NDBO  buoy  measurements.  For  Anita,  the  buoy-derived 

wind  field  represented  conditions  during  a 36-hour  period;  for  Eloise 

the  buoy-derived  wind  field  represented  conditions  during  a 22-hour 

period.  The  radar  data,  however,  in  both  cases  were  collected  in  less 

than  45  minutes.  The  relative  agreement  between  the  buoy-  and  WARF- 

derived  wind-direction  data  was  better  than  5°  for  Eloise  and  20°  for 

Anita.  Comparison  of  the  WARF  and  the  buoy  measurements  of  the  wind 

direction,  coincident  in  both  time  and  space,  agreed  to  within  10°. 

Estimates  of  wind  speed  were  derived  from  a parametric  wave  prediction 
14-17 

model,  which  required  a WARF  measurement  of  wave  height  and  radial 

distance  from  the  center  of  the  storm.  Only  a few  in  situ  wave  height 
and  surface  wind  speed  measurements  were  available  at  the  NDBO  buoys  to 
compare  with  the  WARF  estimates. 

The  WARF  measurements  of  significant  wave  height  and  wind  speed 

made  during  hurricane  Anita  (Appendix  D)  were  compared  to  those  made 

simultaneously  at  NDBO  buoy  EB71  (26.0°N,  93.5°W);  agreement  was  within 

o 1 8 

0.5  m,  0.4  m/s  and  9°  (Appendix  D) . The  methods  used  to  interpret 

rms  wave  height  from  the  sea-echo  Doppler  spectrum  are  accurate  to 
19 

approximately  10%.  The  accuracy  of  the  wind  speed  estimates  are  a 
function  of  the  accuracy  of  the  parametric  wave  prediction  model,  and 


3 


the  WARF  estimates  of  rms  wave  height  and  the  radial  distance.  We  de- 
termined that  the  WARF  measurements  of  significant  wave  height  and 
radial  distance  are  accurate  to  within  0.5  m and  20  km,  respectively. 

For  a significant  wave  height  of  5.5  m and  radial  distances  ranging 

18 

from  30  km  to  300  km,  the  potential  errors  are  less  than  2 m/s.  This 
estimate  is  consistent  with  the  parametric  wave  prediction  model  verifi- 
cation tests  of  Ross  and  Cardone.^ 
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Ill  WARF  SKYWAVE  RADAR 

The  WARF^is  a high-resolution  experimental  high  frequency  (HF) 
skywave  radar  located  in  central  California.  The  radar  is  bistatic 
and  operates  in  the  HF  band  between  6 and  30  MHz.  Ocean  areas  are 
illuminated  by  a 20  kW  swept-f requency  continuous-wave  (SFCW)  signal 
from  a transmitter  site  located  at  Lost  Hills,  California.  The  energy 
reflected  from  the  ocean  surface  is  received  185  km  north  of  the  trans- 
mitting site  at  the  Los  Banos,  California  receiving  site.  The  receiving 
antenna  array  is  2.5-km  long  and  consists  of  a double  linear  array  of 
256  whip  antennas  producing  a nominal  0.5°  azimuthal  beamwidth  at  15  MHz. 
The  signal  propagates  to  and  from  remote  ocean  patches  up  to  3000  km 
from  the  radar  by  means  of  a single  reflection  from  one  or  more  layers 
in  the  ionosphere.  These  layers  are  defined  by  peaks  in  the  electron 
density  height  profile  and  occur  at  different  heights  between  100  and 
500  km. 

The  WARF  coverage  area  is  shown  in  Figure  1 of  Appendix  B.  The 
radar  can  be  directed  either  east  or  west,  and  can  be  electronically 
steered  in  azimuth  +32°  from  boresight  anywhere  within  the  coverage 
area  in  0.25°  increments.  Position  accuracy  is  a function  ot  ionospheric 
height  estimates  where  uncertainties  in  the  midpath  height  result  in  a 
nominal  position  accuracy  of  approximately  20  km.  However,  at  any  one 
location,  the  accuracy  between  consecutive  measurements  in  range  and 
azimuth  is  an  order  of  magnitude  better.  WARF  has  multiple-beam  capa- 
bility. For  hurricane  monitoring,  sea  backscatter  is  simultaneously 
received  at  four  adjacent  ocean  areas  from  four  different  beams  separated 
by  0.25°.  The  size  of  the  ocean  scattering  patch  is  a function  of  the 
radar  beamwidth  and  the  range  cell  separation.  The  size  of  the  minimum 
scattering  patch  is  3 km  in  range  by  15  km  in  azimuth. 


IV  THE  SEA  ECHO  DOPPLER  SPECTRUM 


The  sea  backscatter  received  at  the  WARF  is  coherently  processed  in 

range  and  Doppler  to  produce  a sea-echo  Doppler  spectrum.  The  typical 

sea-echo  Doppler  spectrum  is  characterized  by  two  dominant  first-order 

20 

echoes  surrounded  by  a second-order  continuum.  Crombie  interpreted 
the  first-order  echoes  in  terms  of  simple  Bragg  scattering  that  repre- 
sented a resonant  response  between  radio  waves  of  wavenumber  kQ  and 

ocean  waves  of  wavenumber  k = 2k  . The  radar  measures  the  relative 

o 

power  and  Doppler  of  the  ocean  waves  traveling  radially  toward  or  away 
from  the  radar.  The  power  ratio  of  the  two  first-order  echoes  repre- 
senting the  upwind  and  downwind  conditions,  relative  to  the  radar  look 
direction,  are  indicative  of  the  wave  direction  of  the  waves  with  wave- 
number  k.  Because  k is  usually  large  (k  > 0.5),  it  is  assumed  that  the 
predominant  wind  direction  is  identical  to  the  predominant  direction  of 
these  ocean  waves. 

21  2 2 

Barrick  A»  derived  theoretical  expressions  that  accurately  de- 
scribe the  HF  scattering  process  to  second-order.  For  a specified  di- 
rectional wave  spectrum,  his  model  accurately  predicts  the  sea-echo 
Doppler  spectrum.  Using  Barrick' s second-order  model,  we  have  studied 
the  effects  of  the  wind  direction,  wave  directionality,  and  the  wave 
frequency  spectrum  on  the  Doppler  spectrum.*  The  wave  height  spectrum 
is  derived  from  the  second-order  structure  surrounding  the  first-order 

echoes.  As  wave  energy  increases,  the  amplitude  of  the  second-order 

23 

structure  increases.  Barrick  and  Lipa  have  reviewed  numerous  inver- 
sion techniques  used  to  compute  rms  wave  height,  the  wave  spectrum, 
and  the  directional  wave  spectrum.  Agreement  between  ground-wave 
radar  measurements  and  in  situ  observations  is  within  10%. 


V 


THE  RADAR  MEASUREMENTS 


For  each  coherent  radar  measurement,  we  processed  21  independent 
Doppler  spectra  spaced  contiguously  at  3-km  range  intervals.  These 
spectra  were  obtained  simultaneously  at  each  of  the  four  adjacent  radar 
beams.  The  size  of  the  ocean  patch  monitored  by  the  WARF  was  specified 
by  the  radar  operator.  The  size  was  consistent  with  ionospheric  propa- 
gation conditions  and  the  characteristic  length  scales  of  the  ocean 
wave  and  wind  fields  being  measured.  For  hurricane  sampling,  the  data 
could,  therefore,  be  averaged  over  a large  number  of  adjacent  range  and 
azimuth  cells. 

To  obtain  wind  direction  estimates  with  WARF,  we  used  a coherent 
integration  time  of  12.8  s and  averaged  the  sea-echo  Doppler  spectra 
obtained  from  five  contiguous  range  cells  at  each  azimuth.  For  each 
12.8  s we  obtained  16  wind  direction  estimates.  At  a range  of  3000  km, 
the  typical  distance  of  our  hurricane  measurements,  the  size  of  the 
scattering  patch  is  approximately  15  km  in  range  and  25  km  in  azimuth. 

To  estimate  wave  height  and  wind  speed  with  WARF,  we  used  a coherent 
integration  of  102.4  s and  averaged  the  sea-echo  Doppler  spectra  ob- 
tained from  21  range  cells  and  3 azimuths.  At  a range  of  3000  km,  the 
size  of  the  scattering  patch  is  approximately  60  km  in  range  and  50  km 
in  azimuth.  We  conducted  some  simple  simulations  to  determine  the 
effect  of  the  scattering  patch  size  on  the  radar  measurement  of  wind 
speed  and  direction.  We  only  found  significant  discrepancies  between 
the  simulated  spatial  estimate  and  the  midpoint  value  when  the  scatter- 
ing patch  was  centered  directly  in  the  hurricane  eye.  From  the  region 
of  the  maximum  winds  radially  outward,  we  found  the  relative  agreement 
between  the  space  average  and  the  point  measurement  was  within  approxi- 
mately 10%. 

The  accuracy  of  the  skywave  radar  estimates  of  wave  height  and  wind 
velocity  Is  not  only  dependent  on  the  accuracy  of  the  theoretical 
models,  but  on  the  quality  of  the  measured  sea-echo  Doppler  spectrum. 


The  quality  of  the  recorded  sea  backscatter  depends  on  the  ionospheric 
conditions  over  fairly  short  periods  (tens  of  seconds).  If  the  radio 
waves  propagate  through  a strong,  single,  stable,  coherent  ionospheric 
layer,  than  high-quality  sea  backscatter  similar  to  ground-wave  radar 
is  obtained.  However,  if  the  signals  are  received  at  the  same  time  from 
two  or  more  different  paths  (multipath),  the  backscatter  from  one  ocean 
patch,  generally  shifted  in  Doppler  and  amplitude,  will  be  reflected 
from  a different  part  of  the  ocean  and  a different  part  of  the  iono- 
sphere. The  backscatter  from  the  first  ocean  patch  will  be  contaminated 
by  the  backscatter  from  the  other  ocean  patches.  If  the  ionosphere  is 
al?o  changing  in  time  or  space  during  a coherent  sampling  period, 
smearing  of  the  data  will  occur.  The  ability  to  predict  the  iono- 
spheric conditions  over  these  short  time  periods  would  enhance  the 
quality  of  sea  backscatter  and  reduce  the  sampling  time.  Both  vertical 
and  oblique  incidence  ionospheric  soundings  are  ‘used  to  provide  some 
data  quality  information.  More  time  is  required  to  complete  an  iono- 
spheric sounding  than  is  required  to  take  the  data;  thus,  assessment  of 
the  data  quality  is  difficult  for  rapidly  changing  ionospheres.  There- 
fore, real-time  output  of  the  data  from  the  WARF  site  minicomputer  is 
required  and  is  used  to  verify  the  data  quality. 

The  wind  direction  measurement  is  not  extremely  sensitive  to  iono- 
spheric contamination  caused  by  ionospheric  smearing  and  multipath 
because  only  the  amplitude  of  the  two  strong  first-order  echoes  must  be 
measured,  and  the  coherent  integration  time  is  short  (12.8  s).  Mapping 
the  wind-direction  field  in  a hurricane  can  be  routinely  accomplished 
under  most  ionospheric  conditions  in  10  minutes  or  less.  Once  the 
surface-wind-direction  map  is  made,  the  storm  center  can  be  identified 
and  more  extensive  monitoring  of  the  wind  speed  and  wave  height  anywhere 
within  the  storm  can  be  accomplished. 

The  wave  height  and  wind  speed  estimates,  dependent  on  extracting 
spectral  information  from  the  weak  second-order  echoes  surrounding  the 
first-order  echoes,  are  more  sensitive  to  ionospheric  contamination  than 
the  wind  direction  measurements.  This  contamination  is  the  largest  source 
of  error  in  the  measurements. 
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The  net  effect  of  the  multipath  and  smearing  is  to  increase  the 

calculated  wave  height  and  wind  speed  because  power  is  added  in  an 

unknown  way  to  the  sea-echo  Doppler  spectrum.  Fortunately,  for  the 

large  sea-state  conditions  normally  encountered  during  hurricanes, 

the  power  of  the  second-order  echoes  is  large  and  is  frequently  stronger 

than  the  contamination.  Under  these  conditions,  wave  height  and  wind 

speed  estimates  can  still  be  made.  The  effects  of  contamination  are 

variable  in  time  and  space.  Often  a small  portion  of  the  data  may  not 

be  contamintated  while  the  remainder  of  the  data  will  be  contaminated. 

We  have  developed  several  techniques  to  filter  the  contaminated  sea- 

echo  spectra  out  of  the  sample  to  improve  data  quality.  On-going  work 

at  SRI  and  NOAA  to  develop  better  ways  of  collecting  and  processing 

24 

skywave  data  has  improved  the  skywave  estimates  of  wave  height.  Be- 
cause of  the  longer  integration  times  and  special  data  sampling  and  pro- 
cessing, the  time  required  to  make  an  estimate  of  the  wave  height  and 
wind  speed  is  longer  than  that  necessary  to  estimate  wind  direction. 
However,  the  radar  operator  can  selectively  steer  the  radar  beams  to 
specific  regions  of  the  storm  to  make  the  desired  measurements. 


11 


I 


VI  SUMMARY 

We  have  demonstrated  that  high-resolution  HF  skywave  radar  can  be 
used  to  track  and  to  monitor  remotely  wind  direction  and  speed  and 
wave  height  throughout  all  regions  of  hurricanes  and  tropical  storms. 
The  radar  position  fixes  are  accurate  to  20  km;  the  wind  and  wave 


estimates  are  accurate  to  within  10%.  Skywave  radar  can  be  steered  by 
the  operator  in  real-time  to  specific  regions  of  the  storm  located  over 
3000  km  from  the  radar,  or  the  radar  can  be  used  to  scan  across  the 
entire  storm.  The  skywave  radar  surface  measurements  can  be  combined 
with  the  upper  level  satellite  measurements  to  improve  the  quality  and 
reliability  of  hurricane  and  tropical  storm  observations. 
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REMOTE  MEASUREMENT  OF  THE  POSITION  AND  SURFACE  CIRCULATION 
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ABSTRACT 

The  position  and  surface  wind  direction  field  of  hurricane  Eloise 
was  observed  by  measuring  sea  backscatter  from  the  Gulf  of  Mexico 
through  the  use  of  an  HF  skywave  radar  situated  in  central  California. 

A radar  map  of  the  surface  wind  direction  field  was  compiled  from  the 
direction  of  the  ocean  gravity  waves  8-m  long,  and  compared  to  the  sur- 
face wind  directions  measured  at  National  Data  Buoy  Office  (NDBO)  open 
ocean  moored  buoys  EB10  and  EB04.  Agreement  to  within  10%  was  found. 
The  radar  position  fix  was  determined  from  the  radar-derived  surface 
wind  map  and  showed  agreement  to  within  35  km  of  the  position  estimated 
from  the  official  National  Hurricane  Center  track. 


I INTRODUCTION 

The  position  of  a hurricane  in  open  water  often  can  be  determined 
remotely  by  the  cloud  cover  in  satellite  photographs  and  by 
reconnaissance  aircraft  flights  directly  through  the  center  of  the 
hurricane.  Measurements  of  surface  wind  directions  are  sparse.  We 
discuss  a new  technique  to  estimate  the  surface  wind  directions  and 
position  of  a hurricane  through  the  use  of  high-frequency  (HF)  skywave 
radar  sea  backscatter  measurements.  The  position  and  surface  wind  di- 
rection field  of  hurricane  Eloise  was  determined  remotely  (at  a distance 
of  3000  km)  at  approximately  2030Z  on  22  September  1975  from  a detailed 
description  of  the  ocean  gravity  wave  field. 
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II  THEORY  RELATING  WIND  DIRECTION  TO  RADAR  SPECTRA 


Backscatter  measurements  obtained  from  the  skywave  radar  measure- 
ments are  coherently  range  and  Doppler  processed  to  obtain  a Doppler 
spectrum.  The  sea-echo  Doppler  spectrum  shown  in  Figure  1 is  character- 
ized by  two  dominant  first-order  echoes.  These  strong  echoes  were  first 
explained  by  Crombie  (1955),  and  are  caused  by  the  resonant  interaction 
of  ocean  waves  of  length  L equal  to  one  half  the  electromagnetic  wave  of 
length  X.  Only  ocean  waves  traveling  directly  toward  or  away  from  the 
radar  are  detected  by  this  means.  The  principle  direction  of  ocean 
gravity  waves  (8  to  9 meters  long)  is  inferred  from  the  power  ratio  of 
the  approaching  and  receding  waves.  Because  the  waves  are  short  it  is 
assumed  that  they  are  coincident  with  the  mean  wind  direction.  Tech- 
niques to  calculate  direction  have  been  described  in  the  literature  by 
Long  and  Trizna  (1973),  Ahearn  et  al.  (1974),  Tyler  et  al.  (1974), 

Stewart  and  Barnum  (1975),  and  Barnum  et  al.  (1975).  Stewart  and  Barnum 
(1975)  derived  an  expression  for  the  direction  of  the  ocean  gravity  waves 
based  on  ocean  wave  directional  distribution,  G(0),  of  the  form 


G(0) 


s 

cos 


0_ 

2 


(1) 


suggested  by  Longue t-Higgins , Cartwright,  and  Smith  (1963)  and  Cartwright 
(1963).  The  wave  direction,  6,  relative  to  the  radar  bearing  is 


6 


2 arctan  (r^S) 


(2) 


G(0\ 

where  r = is  the  power  ratio  computed  from  the  first-order 

echoes,  and  s is  a parameter  in  Equation  (1)  which  describes  the 
spreading  envelope  of  the  ocean  gravity  wave  field.  Stewart  and  Barnum 
(1975)  found  that  the  standard  deviation  of  the  measurement  was  ±16  de- 
grees for  open  ocean  conditions.  Equation  (2)  was  used  for  the  data  to 
be  described  herein. 
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DOPPLER  - Hz 


FIGURE  1 


AN  EXAMPLE  OF  A DOPPLER  SPECTRUM  OBTAINED  FROM  PROCESSING 
12.8  s OF  HF  SKYWAVE  RADAR  BACKSCATTER  DATA.  The  first-order 
Bragg  lines  are  the  two  maxima. 
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Ill  BACKSCATTER  SPECTRUM  RECORDING  AND  INTERPRETATION 


The  Wide  Aperture  Research  Facility  (WARF)  HF  skywave  radar  was 
used  to  collect  the  hurricane  data  as  Eloise  moved  northward  through 
the  Gult  of  Mexico.  The  WARF  radar  is  located  in  central  California 
(Figure  2)  approximately  3000  km  west  of  the  hurricane.  The  receiving 
array  is  2.55  km  long  and  forms  a 0.5°  beam  at  15  MHz.  The  area  was 

k 

illuminated  by  a 30-degree  transmitter  beam  centered  on  the  hurricane, 
and  backscatter  energy  on  each  coherent  26-second  radar  measurement  was 
received  at  five  contiguous  one-half  degree  beams  spaced  by  one-quarter 
degree.  The  absolute  position  accuracy  of  each  WARF  dwell  observation 
is  approximately  20  kilometers  in  range  and  20  kilometers  in  azimuth. 

The  relative  accuracy  between  two  measurements  is . significantly  better. 

The  sea  backscatter  was  processed  in  range  and  in  Doppler  to  pro- 
duce a sea-echo  Doppler  spectrum.  Twenty-one  indivudual  Doppler  spectra 
processed  at  21  contiguous  range  cells  were  processed  for  each  12.8  s 
coherent  radar  measurement.  The  Doppler  spectra  from  two  consecutive 
12.8  s periods  were  averaged.  The  Doppler  spectrum  processed  for  each 
range  resolution  cell  was  26  km  in  azimuth  and  3 km  in  range.  An 
example  of  Doppler  spectra  obtained  simultaneously  at  five  different 
azimuths  in  the  vicinity  of  the  eye  of  hurricane  Eloise  is  presented 
in  Figure  3.  These  spectra  are  produced  in  real  time,  and  each  of  the 
five  spectra  represents  a spatial  average  of  63  km  in  range  and  25  km 
in  azimuth.  The  abrupt  change  in  the  first-order  echoes  and,  therefore, 
wind  direction  occurs  over  less  than  ’00  km,  the  typical  diameter  of  a 
hurricane  eye.  At  each  azimuth,  we  averaged  the  Doppler  spectra  from  5 
contiguous  range  cells  (instead  of  21).  Wind  direction  was  estimated 
from  the  Doppler  spectrum  averaged  over  five  range  cells  and  two  time 

periods.  This  procedure  involved  a spatial  average  of  the  wave  direc- 

2 

tion  over  only  390  km  of  ocean  (15  km  in  range  by  26  km  in  azimuth). 
Each  such  measurement  was  centered  at  grid  points  spaced  by  13  km  in 

* Normally,  the  transmitter  beam  is  6°.  Because  one  of  the  transmitters 
was  under  repair  during  the  experiment,  only  the  center  four  elements 
of  the  eighteen  element  transmitting  array  were  used. 
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FIGURE  2 THE  AREA  SCANNED  BY  THE  WIDE  APERTURE  RESEARCH  FACILITY 
(WARF),  SKYWAVE  RADAR  TO  LOCATE  HURRICANE  ELOISE  IN  THE 
GULF  OF  MEXICO  AT  2100Z  ON  SEPTEMBER  22,  1975 
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FIGURE  3 AN  EXAMPLE  OF  THE  DOPPLER  SPECTRA  OBTAINED  FROM  AN  INCOHERENT 
AVERAGE  OF  TWO  COHERENT  12.8  s HF  RADAR  MEASUREMENTS  ACROSS 
THE  EYE  OF  HURRICANE  ELOISE.  Each  Doppler  spectrum  is  an  average  of  the 
Doppler  spectra  obtained  from  21  different  ranges  spaced'  3 km  apart.  Each  Doppler 
spectrum  shown  is  spaced  by  0.25°  in  azimuth.  The  5 Doppler  spectra  were  obtained 
simultaneously  over  an  area  approximately  100  km  in  azimuth  and  60  km  in  range. 
The  amplitude  of  the  first-order  echoes  changes  as  a function  of  wave  direction  from 
approaching  to  receding  directions  in  less  than  100  km. 


azimuth  and  15  km  in  range.  A total  of  20  wind  direction  estimates  are 
obtained  every  26  s.  The  total  time  required  to  complete  a scan  of  the 
entire  hurricane  was  33  min. 

We  used  Equations  (1)  and  (2)  to  derive  the  ocean  gravity  wave 
field  measured  by  the  skywave  radar  for  ocean  wave  lengths  ranging  from 
8 m to  9 m for  three  values  of  s.  For  hurricane  waves,  the  value  of  s 
is  not  well  known.  The  value  of  s is  inversely  proportional  to  wave 
frequency  and  wind  speed  as  shown  by  Ewing  (1969),  Tyler  et  al.  (1974) 
and  others.  Stewart  and  Barnum  (1975)  derived  an  expression  for  s as  a 
function  of  wind  speed  and  wave  number.  Using  their  model  for  hurricane 
wind  speeds  resulted  in  values  of  s less  than  1.0.  This  would  imply  a 
highly  nondirectional  wave  field.  We  calculated  three  wind  direction 
maps  using  s equal  to  1.0,  2.0,  and  4.0. 

The  left-right  ambiguity  in  direction  was  resolved  from  the  general 
circulation  of  a hurricane  in  the  following  manner.  First,  we  assumed 
that  the  directional  properties  of  the  wind  field  and  the  wave  field  for 
the  higher  wave  numbers  were  similar.  Then, we  assumed  that  the  winds 
are  directed  counterclockwise  and  flow  inward,  and  resolved  the  left/ 
right  ambiguity.  Only  in  cases  in  which  the  winds  were  approximately 
radial  to  the  radar  beam  direction  with  no  clear  inflow  angle  discern- 
able  was  the  choice  arbitrary.  In  this  case,  the  possible  direction 
error  was  inconsequential  because  the  difference  in  total  direction  be- 
tween the  left/right  direction  vectors  was  less  than  the  expected  error 
of  the  measurement  (±16  degrees).  The  three  wind  direction  maps  were 
compared  to  the  NWS  surface  analyses  for  1800Z  on  22  September  1975, 
the  two  wind  observations  at  NOAA  buoys  EB10  and  EB04,  and  cloud  pat- 
terns obtained  from  the  NASA  SMS1  Synchronous  Meteorological  Satellite 
photographs.  The  best  agreement  was  found  for  the  wave  pattern  derived 
from  a value  of  s equal  to  1.0.  This  value  of  s was  thus  assumed 
correct  for  all  radar  measures  of  direction  in  the  hurricane. 


27 


IV  HURRICANE  WIND  AND  WAVE  FIELD 


The  directional  ocean  wave  field  derived  from  the  radar  measure- 
ments during  Hurricane  'Eloise  is  shown  in  Figure  4.  During  the  period 
22-23  September  1975,  Hurricane  Eloise  moved  northeasterly  across  the 
Gulf  of  Mexico,  passing  within  56  km  of  NDBO  buoy  EB04  and  within  46  km 
of  buoy  EB10.  We  computed  a wind  direction  field  from  the  buoy-measured 
wind  directions  over  a 22-hr  period  to  compare  to  the  WARF-derived  wind 
directions.  We  plotted  the  buoy-measured  wind  directions  with  respect 
to  one  center.  We  superimposed  the  buoy-derived  wind  field  on  the  WARF- 
derived  wind  map.  The  buoy  data  overlaps  the  WARF  measurements  taken  at 
approximately  2030Z  and  spans  the  period  1000Z  on  22  September  1975 
through  0800Z  on  23  September  1975.  The  agreement  between  the  WARF  and 
buoy  wind  direction  measurements  was  dependent  on  the  assumptions  that 
the  hurricane  surface  wind  circulation  and  storm  motion  remained  uni- 
form during  the  period  when  Eloise  passed  between  the  two  buoys.  How- 
ever, Eloise  stalled  for  a period  of  approximately  5 hours,  so  some 
error  in  both  the  position  and  wind  direction  was  expected.  Nevertheless, 
the  radar  measured  wind  directions  agreed  to  within  10°  overall  with 
the  buoy-derived  spatial  wind  field,  and  they  verified  in  a qualitative 
sense  the  observed  surface  circulation  inferred  for  the  radar  measure- 
ments. The  difference  between  the  wind  directions  measured  at  the  NOAA 
buoys  EB10  and  EB04  (Withee  and  Johnson,  1975)  coincident  in  time  and 
space  with  the  radar-derived  mean  wind  directions  at  those  locations 
was  less  than  10  degrees. 

We  considered  several  possible  sources  of  error  in  the  radar  mea- 
surements. First,  the  full  transmitting  array  was  not  used  during  the 
experiment,  which  resulted  in  a 30  degree  beam  instead  of  the  usual  6 
degree  beam.  Hence  the  combined  radar  antenna  sidelobes  were  higher, 
and  could  have  received  measurable  backscatter  energy  from  other  areas 
of  the  Gulf.  Serious  sidelobe  contamination  was  unlikely  since  the 
sidelobe  gain  of  the  narrow  receive  beam  should  have  been  at  least 
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FIGURE  4 COMPARISON  OF  WARF-MEASURED  SURFACE  CIRCULATION  WITH  NDBO 
DATA  BUOY  AND  SATELLITE  CLOUD  WIND  DIRECTIONS  FOR  HURRICANE 
ELOISE 
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20  dB  down.  Backlobe  reception  of  scatter  from  the  Pacific  Ocean  was 
also  measured  by  reversing  the  transmitter  steer,  and  it  was  found  to 
be  negligible  during  the  experiment.  Second,  ionospheric  multipath, 
the  simultaneous  reception  of  signals  from  two  different  patches  of  the 
ocean  could  contaminate  the  sea-echo  Doppler  spectrum.  Especially, 
near  the  storm  center  where  the  wind  direction  changes  significantly 
over  small  distances,  multipath  could  cause  errors  in  computing  the 
power  ratio  of  the  first-order  echoes.  By  carefully  managing  the 
propagation  conditions  with  real-time  vertical-  and  oblique-incidence 
soundings  and  visual  inspection  of  the  data  for  quality,  we  minimized 
the  effects  of  multipath  for  Eloise. 


V POSITION  OF  HURRICANE  ELOISE 

The  position  of  the  hurricane  was  determined  from  the  direction 
field  derived  for  each  of  the  three  values  of  s.  The  three  positions 
were  internally  consistent  to  within  10  km.  The  position  of  Eloise, 
shown  in  Figure  4 for  s = 1.0  is  within  35  km  in  range  and  aligned  in 
azimuth  with  the  position  determined  from  the  NHC  track.  The  accuracy 
of  determining  the  position  of  the  hurricane  is  primarily  dependent  on 
the  position  accuracy  of  the  WARF  and  the  interpretation  of  the  wave 
field  near  the  eye.  The  position  accuracy  of  the  WARF  is  dependent  on 
the  estimate  of  the  virtual  height  of  the  ionosphere  at  the  time  of  the 
measurement.  Ionospheric  soundings  or  a reference  echo  received  from  a 
known  location  can  improve  the  estimate  of  virtual  height.  The  accuracy 
of  the  WARF  using  ionospheric  soundings  is  20  km.  When  we  used  a high- 
frequency  repeater  for  measurements  along  the  Gulf  coast,  the  accuracy 
of  the  observations  could  be  increased  to  5 km  or  better.  The  left/ 
right  ambiguity  in  the  direction  of  the  wave  field  near  the  eye  of  the 
hurricane  is  difficult  to  resolve  without  prior  knowledge  of  the  po- 
sition of  the  eye.  This  center  position  is  determined  from  scans 
similar  to  Figure  3.  The  left/right  ambiguity  is  then  resolved.  For 
the  Eloise  data,  this  ambiguity  could  cause  a 10  to  30  km  error  in 
our  position  of  the  hurricane. 
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V 


CONCLUSION 


The  position  and  surface  circulation  of  the  wind  field  was  deter- 
mined remotely  by  HF  skywave  radar  for  Hurricane  Elolse  from  backscatter 
measurements  of  the  direction  of  the  short  ocean  gravity  waves.  The 
position  of  the  hurricane  agreed  to  within  35  km  of  the  reported  posi- 
tion of  the  hurricane.  The  circulation  agreed  with  that  estimated 
using  NWS  surface  charts,  SMS1  satellite  cloud  photographs,  and  two  NOAA 
buoy  measurements.  The  ability  to  infer  detailed  wind  direction  in  the 
vicinity  of  the  eye  of  the  hurricane  and  in  the  region  of  maximum  winds 
was  demonstrated  and  suggests  a potential  for  measuring  the  waves, 
currents,  and  winds  at  the  surface  level  in  this  region. 
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TRACKING  HURRICANE  ANITA  BY  HF  SKYWAVE  RADAR 

Joseph  W.  Maresca,  Jr.  and  Christopher  T.  Carlson 
SRI  International,  Menlo  Park,  CA 


Hurricane  Anita  was  tracked  across  the  Gulf  of  Mexico  from  29  August 
1977  to  2 September  1977  by  using  an  HF  skywave  radar.  The  radar 
position  estimates  were  made  from  California,  3 to  5 times  per  day,  at 
the  SRI-operated  Wide  Aperture  Research  Facility  (WARF) . The  storm 
center  was  located  using  the  WARF-derived  ocean  surface  wind  direction 
maps.  Seventeen  independent  wind  maps  were  used  to  develop  the  WARF- 
derived  track.  The  WARF-derived  positions  were  compared  to  coincident 
temporal  positions  on  the  official  track  produced  by  the  National 
Hurricane  Center  (NHC) , and  the  relative  agreement  was  19  km. 


'’’RACKING  HURRICANE  ANITA  BY 
HF  SKYWAVE  RADAR 

Joseph  W.  Maresca,  Jr.  and  Christopher  T.  Carlson 

SRI  International 
Menlo  Park,  California  94025 


INTRODUCTION 

This  note  describes  a new  method  for  tracking  hurricanes  using  an 
HF  skywave  radar.  Hurricane  Anita,  the  first  Atlantic  hurricane  of  the 
1977  season,  formed  as  a tropical  depression  in  the  Gulf  of  Mexico  at 
about  1200Z  on  29  August  1977.  Anita  developed  into  a tropical  storm 
at  approximately  0600Z  on  30  August  1977  and  intensified  into  a hurri- 
cane about  12  hours  later.  As  the  storm  moved  west  across  the  Gulf, 
maximum  winds  in  excess  of  75  m/s  were  recorded  prior  to  landfall 
approximately  48  km  south  of  Brownsville,  Texas. 

The  skywave  radar  measurements  of  Anita  were  made  at  the  California- 
based  Wide  Aperture  Research  Facility  (WARF)  [SRI,  1977],  an  SRI- 
operated  high-resolution  HF  skywave  radar.  Seventeen  independent  posi- 
tion estimates  were  computed  for  Hurricane  Anita  from  WARF  at  ranges 
between  2400  km  and  3200  km,  beginning  on  29  August  1977  and  continuing 
through  landfall  on  2 September  1977.  The  hurricane  position  estimates 
were  derived  from  skywave  radar  maps  of  the  areal  distribution  of  the 
surface  wind  directions  generated  using  techniques  similar  to  those  used 
for  Hurricanes  Caroline  [Maresca  and  Barnum,  1975],  Eloise  [Maresca, 

1976  and  Maresca  and  Barnum,  1978],  Kate  [Maresca  and  Carlson,  1977], 
and  Babe  [Maresca  and  Carlson,  1977,  1978].  The  Anita  maps  were  updated 
3 to  5 times  per  day  at  arbitrary  times  including  nighttime  periods. 

The  WARF  radar  position  estimates  cluster  around  the  corresponding 
interpolated  temporal  positions  on  the  National  Hurricane  Center  (NHC) 
official  smooth  track.  The  NHC  track  was  obtained  from  satellite 
Imagery,  reconnaissance  aircraft,  shore-based  microwave  Doppler  radar, 
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and  open  ocean  moored  buoys.  The  mean  deviation  of  the  WARF  positions 
from  corresponding  temporal  positions  along  the  NHC  track  is  19  km. 

WARF  skywave  radar  position  estimates  for  Hurricanes  Caroline,  Eloise, 
Kate,  and  Babe  show  agreement  with  NHC  positions  ranging  from  25  to  40 
km.  The  number  of  WARF  position  estimates  for  each  of  these  hurricanes 
was  limited  to  only  one  or  two.  We  attribute  the  improved  agreement 
between  the  WARF-  and  NHC-derived  tracks  to  the  increased  number  of 
WARF  wind  maps  that  showed  the  evolving  wind  circulation  with  time. 

The  accuracy  of  the  skywave  radar  position  estimates  is  similar  to  the 
accuracy  of  the  satellite  position  estimates  [Sheets  and  Grieman,  1975]. 
One  can  more  reliably  track  hurricanes  remotely  by  combining  both  sky- 
wave  radar  and  satellite  cloud  measurements,  thus  reducing  the  frequency 
of  aircraft  reconnaissance  flights  required  principally  for  identifica- 
tion of  the  hurricane  center. 

Skywave  measurements  of  significant  wave  height,  ocean  wave  spectra 
[Maresca  and  Carlson,  1978;  Maresca,  1978;  Maresca  and  Georges,  1978], 
and  the  mean  surface  wind  velocity  [Maresca  and  Carlson,  1977,  1978] 
were  also  made  for  Hurricanes  Anita  and  Babe  and  were  compared  to  mea- 
surements made  at  NOAA  Data  Buoy  Office  (NDBO)  data  buoy  EB71.  Agree- 
ment between  the  radar-  and  buoy-measured  wind  speeds  was  within  +2  m/s 
(10%),  and  agreement  between  the  wind  direction  measurements  was  within 
10°.  The  radar  estimates  of  significant  wave  height  were  within  0.5  m 
(10%)  of  the  buoy  measurements. 


WARF  SKYWAVE  RADAR 

WARF  is  an  SRI-operated  experimental  bistatic  HF  skywave  radar 
located  in  Central  California.  The  radar  is  operated  in  the  high- 
frequency  (HF)  band  between  6 and  30  MHz.  A 20-kW  swept-frequency  con- 
tinuous-wave (SFCW)  signal  is  transmitted  from  the  transmitter  site  at 
Lost  Hills,  California,  and  after  round-trip  ionospheric  propagation, 
is  received  at  a site  185  km  to  the  north  at  Los  Banos,  California. 

The  receiving  array  is  2.5  km  long  and  forms  a beam  of  0.5°  at  15  MHz. 
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The  WARF  radar  can  be  electronically  steered  in  0.25°  increments 
anywhere  within  the  coverage  area  shown  in  Figure  1.  The  radio  waves 
propagate  over  long  distances  to  an  ocean  scattering  patch  via  one  or 
more  ionospheric  "reflections".  Remote  ocean  patches  illuminated  by 
only  one  ionospheric  reflection  are  typically  within  the  minimum  radar 
range  of  1000  km  and  the  maximum  range  for  one-hop  propagation  of  more 
than  3000  km.  The  nominal  absolute  position  accuracy  is  about  20  km 
without  any  surface  position  references,  but  at  any  one  location  the 
accuracy  between  consecutive  measurements  in  range  and  azimuth  is  an 
order  of  magnitude  better. 

The  size  of  the  ocean  patch  monitored  by  the  WARF  radar  is  a func- 
tion of  the  sampling  parameters  specified  by  the  radar  operator.  The 
operator  considers  the  characteristic  length  scales  of  the  ocean  wave 
features  before  selecting  the  size  of  the  ocean  patch  illuminated  by 
the  radar.  For  example,  in  a hurricane  where  sea  conditions  may  change 
rapidly  over  short  distances,  the  monitored  ocean  patch  would  be  smaller 
than  the  one  monitored  for  large  ocean  storms.  We  usually  sample  21 
contiguous  range  cells  spaced  at  3 km.  The  azimuthal  or  cross-range 
resolution  is  a function  of  the  radar  beamwidth.  For  a one-hop  range 
of  3000  km,  the  0.5°  beam  would  produce  a cross-range  distance  of  about 
25  km.  We  have  found  from  simulation  experiments  [Maresca  and  Carlson, 
1977]  that  a scattering  patch  15  km  in  range  by  25  km  in  azimuth  is 
sufficiently  small  to  accurately  determine  the  wind  and  wave  parameters 
of  the  storm. 


RADAR-MEASURED  WIND-DIRECTION  MAPS 

The  HF  skywave  radar-measured  surface  wind  directions  are  derived 
from  the  predominant  direction  of  decameter  ocean  gravity  waves  about 
10  m long,  which  are  assumed  to  be  tightly  coupled  to  the  wind  over  a 
period  of  tens  of  minutes.  Measurements  of  the  directional  wave  spectrum 
are  generally  sparse,  but  available  measurements  [Longuet-Higgins  et 
al.,  1963;  Ewing,  1969;  Mitsuyasu  et  al.,  1975]  support  the  computation 
of  wind  direction  from  the  wave  direction.  The  WARF  radar  measures  the 
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FIGURE  1 COVERAGE  AREA  OF  THE  WARF  HF  SKYWAVE  RADAR.  All  Anita 
measurements  were  made  west  of  88°W  in  the  Gulf  of  Mexico. 
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relative  power  of  the  approaching  and  receding  ocean  waves  that  satisfy 

the  resonance  condition  with  the  radio  waves  k = 2k  , where  k is  the 

o o 

radio  wave  number  and  k the  ocean  wave  number.  An  example  of  the  radar- 
measured  sea-echo  Doppler  spectrum  is  shown  in  Figure  2.  This  figure 
shows  approximately  13  dB  more  energy  in  the  waves  receding  from  the 
radar  than  in  those  approaching  the  radar.  If  the  directional  distri- 
bution of  the  ocean  waves  were  known  precisely,  the  predominant  wind 
direction  could  be  determined  directly  from  the  radar  measurement.  It 
is  possible  to  compute  the  directional  distribution  directly  from  the 
radar  data  using  Lipa's  [ 1977a, b]  inversion  model.  However,  this  method 
is  computationally  difficult,  and  requires  exceptionally  high-quality 
data  produced  from  coherent  radar  dwells  of  100  s or  more.  To  simplify 
our  measurement  of  wind  direction,  we  have  used  several  models  to 
describe  the  directional  wave  distribution.  We  used  a cosine  model 
proposed  by  Longue t-Higgins  et  al.  [1963]: 

G(0)  - cos8  — (1) 


where  G(0)  is  the  directional  distribution,  0 is  the  angle  between  the 
radar  line-of-sight  and  wind  (wave)  direction,  and  s is  a parameter 
describing  the  spread  of  the  waves.  The  radar-measured  power  ratio 
between  the  approaching  and  receding  ocean  waves  is 


G(0) 

G(0-hr) 


(2) 


Solving  for  0 and  substituting  Eq.  (1)  into  Eq.  (2)  gives 

0-2  arctan  r1/s  (3) 


Several  models  have  been  proposed  to  calculate  0 from  the  radar  data 
[Long  and  Trizna,  1973;  Ahearn  et  al.,  1974;  Tyler  et  al.,  1974; 

Stewart  and  Barnum,  1975;  Barnum  et  al.,  1977;  Maresca  and  Carlson, 
1977].  The  accuracy  of  the  0 estimate  is  dependent  on  the  validity  of 
the  directional  model  and  the  value  of  s.  For  nonhurricane  winds, 
estimates  of  the  mean  radar-derived  wind  direction  agree  to  within  ±16° 
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FIGURE  2 AVERAGE  SEA-ECHO  DOPPLER  SPECTRUM  RECORDED  WITHIN  35  km 

OF  THE  CENTER  OF  HURRICANE  ANITA  AT  2343Z  ON  AUGUST  31,  1977 
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of  the  wind  direction  measured  on  research  vessels  and  ships  of  oppor- 
tunity [Stewart  and  Barnum,  1975],  We  have  used  the  methods  derived  by 
Stewart  and  Barnum  [1975]  and  Mitsuyasu  et  al.  [1975]  to  calculate  s 
and  0 Both  methods  are  based  on  the  wind-speed  to  wave-speed  scaling 
where  G(0)  is  broad  at  high  frequencies  and  narrow  at  low  frequencies. 
For  hurricane  winds,  both  methods  predict  values  of  s considerably  less 
than  1.0  which  results  in  a directional  spread  which  is  too  broad. 
Comparison  of  the  radar  measurements  with  coincident  offshore  NDBO  data 
buoy  measurements  suggests  that  appropriate  values  of  s range  from  1.0 
to  2.0.  Regier  [1975]  has  shown  for  a fixed-wave  frequency  that  the  di- 
rectional distribution  narrows  as  wind  speed  increases.  This  supports 
the  larger  values  of  s found  in  our  measurements. 

The  value  of  s used  to  compute  0 is  based  on  our  previous  experi- 
ments [Maresca  and  Barnum,  1975;  Maresca  and  Carlson,  1977]  and  on  in 
situ  spot  measurements  of  the  wind  direction.  We  expect  the  value  of  s 
to  vary  across  the  hurricane,  but  no  attempt  has  been  made  to  account 
for  this  effect.  Nevertheless,  application  of  this  directional  model 
does  result  in  wind  direction  maps  describing  the  hurricane  surface- 
wind  directions  which  agree  to  within  5 to  10°  with  available  in  situ 
measurements  [Maresca  and  Carlson,  1977,  1978),  and  to  within  10  to  20° 
with  buoy-derived  wind  direction  fields  produced  from  data  recorded 
over  20-hour  periods.  The  solution  of  Eq.  (3)  for  0 results  in  a left- 
right  wind  direction  ambiguity  relative  to  the  radar  line  of  sight. 

This  ambiguity  can  be  resolved  from  well-known  cyclonic  patterns  of  a 
tropical  storm. 

SKYWAVE  RADAR-DERIVED  TRACK  OF  HURRICANE  ANITA 

Twenty-one  independent  position  estimates  were  made  for  Hurricane 
Anita  over  a 5-day  period.  The  first  4 radar  position  estimates  were 
not  used  in  the  WARF  track  because  the  storm  had  two  centers  during 
this  early  period  as  we  observed  in  both  the  satellite  cloud  photo- 
graphs and  WARF  wind  maps.  For  the  remaining  17  estimates,  we  updated 


the  position  estimates  3 to  5 times  per  24-hour  period.  The  radar 


position  estimates  are  based  solely  on  the  surface-wind  circulation  maps. 
Examples  of  the  radar-derived  surface  wind  circulation  maps  made  over  4 
days  are  shown  in  Figure  3.  In  Figure  4 a surface  wind  direction  field 
derived  from  data  recorded  by  NDBO  buoy  EB71  is  compared  to  one  of  the 
WARF  maps.  These  buoy-measured  wind  directions  were  recorded  at  2-hour 
intervals  during  the  period  ±18  hours  of  Anita's  passing  EB71.  We  used  a 
time-space  conversion  to  compute  the  buoy-derived  wind  field  by  assuming 
uniform  wind  direction  and  lateral  storm  motion  during  this  period.  We 
compared  the  buoy-derived  wind  directions  to  the  WARF-derived  wind  direc- 
tions; agreement  was  within  19°.  Agreement  between  the  WARF-derived  wind 
direction  estimate  coincident  in  time  and  space  with  the  buoy  wind 
direction  estimate  was  1°.  Each  map  took  between  30  and  60  minutes  to 
complete  although  the  measurements  near  the  center  used  for  positioning 
took  less  than  10  minutes.  A single  radar  measurement  of  wind  direction 
producing  16  direction  arrows  requires  an  integration  time  of  only 
12.8  s;  therefore,  the  entire  wind  map  could  be  produced  in  a 10-min 
period.  Our  maps  took  longer  to  complete  because  we  took  60-s  radar 
dwells  to  allow  postexperiment  analysis  of  wave  height  and  wind  speed. 

Using  maps  similar  to  the  ones  in  Figure  3,  we  developed  a smooth 
skywave  radar  track  for  Anita.  The  WARF  position  estimates  were  com- 
pared to  the  interpolated  positions  estimated  on  the  official  smooth 
track  of  the  NHC  shown  in  Figure  5 . The  NHC  track  was  derived  from  an 
integration  of  individual  position  estimates  obtained  from  satellite 
photographs,  periodic  aircraft  reconnaissance,  shore-based  microwave 
radar  and  NDBO  buoys.  The  mean  displacement  between  the  WARF  position 
estimates  and  the  interpolated  temporal  position  estimates  along  the 
NHC  smooth  track  is  ±19  km. 


MEASUREMENT  ERRORS 

The  precision  of  the  radar-tracking  measurements  is  dependent  on 
the  absolute  position  accuracy  of  the  WARF  radar,  the  accuracies  assoc- 
iated with  techniques  used  to  produce  the  surface  wind  maps  by  radar, 
and  the  uncertainties  in  determining  a storm  center  based  on  the  surface- 
wind  circulation  field. 
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FIGURE  3 WARF-DERIVED  SURFACE  WIND  MAPS  FOR  ANITA  MADE  AT  (a)  1700Z 
ON  AUGUST  30,  1977,  (b)  2140Z  ON  AUGUST  31,  1977,  (c)  1934Z  ON 
SEPTEMBER  1,  1977,  AND  (d)  0456Z  ON  SEPTEMBER  2,  1977 
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FIGURE  3 (Continued) 
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FIGURE  4 COMPARISON  OF  WARF-DERIVED  SURFACE  WIND  DIRECTION  MAP 

MADE  FOR  ANITA  BETWEEN  2126Z  AND  2156Z  ON  AUGUST  31,  1977 
AND  EB71-DERIVED  SURFACE  WIND  DIRECTION  MAP  MADE  FOR 
ANITA  BETWEEN  0600Z  ON  AUGUST  31,  1977  AND  1800Z  ON 
SEPTEMBER  1,  1977 
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FIGURE  5 COMPARISON  OF  THE  WARF-DERIVED  POSITIONS  WITH  THE  NHC-DERIVED 
TRACK  OF  HURRICANE  ANITA 


47 


A. 


Absolute  Position  Accuracy  of  WARF 


The  absolute  position  accuracy  of  the  WARF  radar  depends  in  part 
on  the  precision  of  midpath  ionospheric  height  estimates.  However, 
direct  measurements  of  midpath  ionospheric  height  are  generally  imprac- 
tical to  make.  We  use  overhead  vertical-incidence  soundings  to  esti- 
mate ionospheric  height.  In  the  absence  of  a reference  echo,  we  can 
obtain  an  absolute  range  accuracy  of  approximately  20  km  using  only  the 
vertical-incidence  ionospheric  height  measurements.  The  ionospheric 
height  estimates  based  on  the  sounding  measurements  are  accurate  to 
within  5 km  to  10  km.  At  a range  of  3000  km,  a 10-km  height  error  for 
propagation  at  a 300-km  F-layer  elevation  produces  a range  error  of 
±6  km.  Ten  years  of  operating  experience  at  WARF  suggests  that  a 20  km 
accuracy  is  typical  unless  ionospheric  propagation  conditions  are  com- 
plex. Estimation  of  ionospheric  height  thus  represents  the  largest 
uncertainty  in  tracking  hurricanes  by  HF  skywave  radar. 

Ionospheric  height  estimates  can  be  improved  significantly  if 
reference  echoes  produced  by  land,  oil  platforms,  or  other  fixed 
structures  are  identified  in  the  radar  data.  These  references  provide 
a natural  means  of  estimating  the  midpath  ionospheric  height. 

In  addition  to  the  range  errors  discussed  above,  errors  in  azimuth 
may  be  caused  by  instrument  imprecision,  coning  errors,  and  ionospheric 
tilts.  We  assume  a mean  absolute  azimuth  error  of  0.5°,  resulting  in  a 
potential  azimuthal  error  of  approximately  25  km.  This  type  of  error 
is  difficult  to  recognize  unless  a reference  echo  is  measured  simultan- 
eously by  the  radar. 

B.  Accuracy  of  Radar-Derived  Surface  Wind  Maps 

The  location  of  the  storm  center  by  radar  is  not  strongly  influenced 
by  the  accuracy  of  the  directional  model.  Reasonable  estimates  of  s, 
and  subsequently  6,  are  sufficient  to  produce  an  internally  consistent 
map  from  which  to  locate  the  storm  center.  The  primary  sources  of  error 
in  constructing  a radar-derived  surface-wind  distribution  map  are  the 
directional  ambiguity  inherent  in  Eq.  (3),  the  influence  of  ocean  scat- 
tering patch  dimensions,  and  the  density  of  the  radar  measurements  near 
the  storm  center. 
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The  method  of  resolving  the  left/right  ambiguity  is  described  in 
Maresca  and  Carlson  [1977].  Near  the  center,  the  effect  of  the  ambi- 
buity  in  azimuth  is  negligible  because  the  power  ratio  of  the  approach- 
ing to  receding  waves  changes  sign.  The  ambiguity  is  more  difficult  to 
resolve  in  range  because  the  wind  direction  is  approximately  perpendic- 
ular to  the  radar  beam.  An  iterative  technique  is  used  to  select  a 
center  beginning  in  the  outer  regions  of  the  storm  where  the  ambiguity 
is  easily  resolved.  Nevertheless,  the  final  selection  of  a center  is 
somewhat  subjective.  After  several  maps  are  compiled,  the  pattern  that 
develops  in  the  surface  wind  circulation  simplifies  the  selection.  The 
error  caused  by  the  ambiguity  is  difficult  to  isolate  and  quantify. 
Further  comparisons  of  surface  wind-direction  measurements  made  simul- 
taneously by  WARF  and  by  other  sensors  will  help  quantify  the  overall 
accuracy  of  the  radar-derived  wind  maps. 

The  dimensions  of  the  ocean  scattering  patch  used  in  these  mea- 
surements was  15  km  in  range  by  25  km  in  azimuth.  The  effect  of  the 
size  of  the  scattering  patch  on  the  wind  speed  and  direction  estimates 
throughout  different  regions  of  the  hurricane  was  investigated  in 
Maresca  and  Carlson  [1977].  The  simulations  suggested  that  only  in  the 
eye  region  (±20  km  of  the  center)  where  a 180°  shift  in  wind  direction 
is  possible  over  the  scattering  patch  does  the  potential  for  large 
errors  exist.  Outside  the  eye,  from  the  region  of  maximum  winds  out- 
ward, the  errors  are  minimal. 

The  measurement  density  or  spatial  coverage  of  the  radar  affects 
the  accuracy  of  the  hurricane  position  estimate.  The  measurements 
made  during  Caroline  and  Kate  were  not  as  dense  as  those  made  for  Eloise, 
Anita,  and  Babe.  As  a result,  the  storm  centers  for  Caroline  and  Kate 
were  less  clearly  defined  than  for  the  other  hurricanes.  The  reliability 
of  the  position  estimate  improved  with  an  increase  in  the  number  and 
density  of  the  measurements.  Experience  dictates  the  importance  of 
completely  mapping  the  entire  wind  pattern  near  the  eye  for  reliable 
position  estimates  based  on  surface  wind  circulation. 

C.  Accuracy  of  Hurricane  Position  Based  on  Surface  Wind  Circulation 

The  radar-derived  surface-wind-direction  maps  are  internally  con- 
sistent with  respect  to  absolute  position.  Comparison  of  corresponding 
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NHC-  and  WARF-derived  estimates  of  hurricane  position  for  the  five 
hurricanes  shows  agreement  ranging  from  5 km  to  50  km,  based  upon 
surface-wind  field  measurements.  Hurricanes  are  highly  variable  in 
time  and  space.  It  is  difficult  to  accurately  model  a "universal" 
hurricane  from  which  the  model  surface-wind  direction  would  provide  an 
accurate  position  estimate  for  every  hurricane.  The  stoms  are  fre- 
quently asymmetric  and  must  be  analyzed  individually,  according  to 
their  specific  characteristics.  An  extensive  analysis  of  the  radar, 
buoy,  and  satellite  wind  fields  has  not  yet  been  undertaken  to  find 
better  ways  of  interpreting  the  radar  data.  A more  extensive  compari- 
son of  individual  radar-  and  buoy-measured  wind  direction  values  will 
lead  to  an  improved  empirical  model. 


SUMMARY 

The  tracking  experiments  conducted  to  date  are  limited  to  5 
different  hurricanes.  Only  the  Hurricane  Anita  data  set  was  sufficient 
to  develop  a smooth  radar-derived  track  spanning  several  days  to  compare 
to  the  NHC  track.  Agreement  between  the  NHC  Best  Track  and  the  WARF- 
measured  track  was  within  ±19  km.  The  WARF  radar  can  be  used  indepen- 
dently or  in  conjunction  with  other  sensors  to  estimate  hurricane 
position. 

The  radar  can  reduce  the  chances  of  large  errors  caused  by  in- 
correctly identifying  the  hurricane  eye  in  the  satellite  cloud  photo- 
graphs, can  confirm  that  the  cloud-level  position  estimate  is  represen- 
tative of  the  surface  position,  can  provide  data  when  the  cirrus  clouds 
completely  obscure  the  storm  center  to  photographic  analysis,  and  can 
provide  additional  surface  analysis  during  the  early  development  stages 
when  the  center  may  be  diffuse  or  when  several  centers  may  be  simul- 
taneously observed.  This  additional  data  should  be  especially  helpful 
for  analyzing  storms  of  weak  to  moderate  strength.  The  radar  position 
estimates  can  be  made  during  both  day  and  night  periods. 
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HF  SKYWAVE  RADAR  TRACK  OF  TROPICAL  STORM  DEBRA 


J.  W.  Maresca,  Jr.  and  C.  T.  Carlson 
SRI  International,  Menlo  Park,  CA 


Tropical  storm  Debra  was  tracked  over  a two  day  period  using  the 
Wide  Aperture  Research  Facility  (WARF)  HF  skywave  radar.  Tracking 
began  at  1800Z  on  27  August  1978  while  Debra  was  only  a depression  and 
continued  through  landfall  at  about  0000Z  on  29  August  1978.  Seven 
position  estimates  were  made  from  radar-derived  surface  wind  direction 
maps.  At  about  1800Z  on  28  August  1978  Debra  intensified  to  a tropical 
storm,  and  the  WARF  radar  position  estimates  were  compared  to  recon- 
naissance aircraft,  the  Galveston  microwave  Doppler  radar,  and  satellite 
position  estimates  for  two  different  time  periods.  Agreement  was 
25±15  km,  22+10  km  and  63±10  km,  respectively. 


HF  SKYWAVE  RADAR  TRACK  OF  TROPICAL  STORM  DEBRA 


J.  W.  Maresca,  Jr.  and  C.  T.  Carlson 

SRI  International 
Menlo  Park,  California  94025 


I INTRODUCTION 

The  surface  wind  and  wave  fields  of  hurricanes  can  be  monitored  re- 
motely using  a high-resolution  HF  skywave  radar.  Propagation  of  the 
radio  waves  via  the  ionosphere  permits  propagation  to  distances  greater 
than  3000  km.  We  conducted  experiments  at  the  SRI-operated  Wide  Aperture 
Research  Facility  (WARF)  HF  skywave  radar  that  demonstrated  the  feasi- 
bility of  tracking  tropical  storms  and  mapping  the  intensity  of  the 
winds  and  waves  in  all  regions  of  the  storm. * ^ The  surface  wind  di- 
rection maps  can  be  routinely  produced.  A track  can  be  developed  by 
locating  the  center  of  the  storm  from  each  wind  direction  map.  A track 
was  developed  for  hurricane  Anita^  from  17  wind-direction  maps.  The 
WARF  radar-determined  positions  were  compared  to  interpolated  positions 
along  the  official  track  developed  by  the  National  Hurricane  Center 
(NHC) , and  the  relative  agreement  was  within  19  km. 

Between  1800Z  on  27  August  1978  and  0100Z  on  29  August  1978,  the 
WARF  skywave  radar  was  used  to  track  tropical  storm  Debra.  Debra  in- 
tensified from  a weak  tropical  depression  to  a tropical  storm  during 
this  time.  The  purpose  of  this  paper  is  to  show  the  results  of  long- 
range  radar  tracking  of  a weak  storm  that  has  a large  and  poorly- 
defined  center.  All  of  the  Debra  measurements  were  made  at  distances 
greater  than  2800  km  from  the  radar  via  the  F layer.  The  radar  position 
estimates  were  computed  from  seven  WARF-derived  surface  wind  direction 
maps,  and  were  compared  to  satellite,  reconnaissance  aircraft,  and  the 
Galveston  shore-based  radar  position  fixes.  All  position  fixes  were 
compared  to  the  aircraft  position  fixes  at  approximately  2100Z  and 


2400Z  on  28  August  1978.  The  mean  differences  between  the  aircraft  po- 
sitions and  the  WARF-radar,  satellite,  and  Galveston  radar  were  25±15  km 
44±6  km,  and  7+3  km,  respectively. 


II  THE  MEASUREMENTS 

A storm  center  is  determined  from  the  WARF-maps  of  the  surface  wind 
direction  field.  The  WARF  wind  direction  estimates  are  derived  from 
the  predominant  direction  of  ocean  waves  over  a rectangular  area  approxi- 
mately 15  km  x 25  km.  The  sea  backscatter  is  range  and  Doppler  pro- 
cessed for  coherent  integration  times  of  12.8  s to  produce  a sea-echo 
Doppler  spectrum.  Wind  direction  is  estimated  from  the  ratio  of  the 
upwind  to  downwind  first-order  echoes  of  the  Doppler  spectrum.  A sur- 
face wind  map  is  produced  from  500  to  1000  individual  wind  direction 
estimates.  The  time  required  to  map  the  wind  direction  field  for 
Debra  was  approximately  45  min.  However,  this  time  can  be  reduced  to 
10  min  under  operational  situations,  because  the  WARF  usually  requires 
less  than  10  min  to  collect  and  analyze  all  the  data  required  to  iden- 
tify a storm  center. 

The  accuracy  of  locating  the  storm  center  from  the  WARF-derived 
wind  maps  is  dependent  on  the  absolute  position  accuracy  of  the  radar, 
the  accuracy  of  the  wind  direction  maps,  and  the  uncertainties  in 
determining  the  storm  center  from  the  surface  wind  circulation.  The 
location  of  the  wind  map  is  dependent  on,  among  other  things,  the 
height  and  tilt  of  the  ionosphere.  Estimates  of  the  ionospheric 
height  can  be  made  from  overhead  vertical-incidence  soundings.  In  some 
instances,  these  soundings  are  not  representative  of  the  actual  propa- 
gation conditions  at  midpath.  Uncertainties  in  properly  estimating  the 
midpath  ionospheric  height  can  result  in  errors  of  ±20  km  or  more.  We 
improved  the  accuracy  of  our  estimates  of  the  ionospheric  height  for 
Debra  over  those  made  for  hurricane  Anita  by  recording  backscatter  along 
the  coast  both  before  and  after  the  data  for  each  wind  map  was  collected. 
The  ionospheric  height  was  calculated  from  the  known  location  of  the 
land  and  this  height  was  used  tc  locate  the  wind  maps. 
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Ionospheric  motions  and  irregularities  and  multiple-path  propaga- 
tion may  contaminate  the  radar  measured  sea-echo  Doppler  spectrum  and 
may  require  more  radar  measurements  to  obtain  high  quality  data  at  each 
measurement  location.  These  additional  radar  measurements  increase  the 
time  required  to  produce  a wind  map.  Severe  contamination  may  prevent 
the  collection  of  data  at  specific  locations.  The  wind  direction  mea- 
surement, however,  is  not  extremely  sensitive  to  this  contamination, 
and  more  than  95%  of  the  data  recorded  for  Debra  from  sequential  radar 
scans  were  used  in  this  analysis.  Moreover,  the  ionosphere  was  severely 
disturbed  dnring  most  of  this  tracking  experiment.  In  addition,  a 
strong  sporadic-E  layer  prevented  propagation  to  the  Gulf  of  Mexico 
between  1730Z  and  2030Z  on  28  August  1978.  This  was  the  first  time  in 
all  six  hurricanes  that  were  tracked  with  WARF  since  1975  that  the 
propagation  conditions  prevented  data  collection.  Nonetheless,  the 
WARF  measurements  of  Debra  were  taken  during  both  day  and  night,  and 
wind  maps  were  made  every  3 to  5 hours  during  our  experiment. 


Ill  WARF-DERIVED  TRACKS 

Four  of  the  seven  WARF  wind  maps  made  for  Debra  are  shown  in  Figure 
1.  The  maps  are  representative  of  the  development  of  the  storm  from  a 
weak  depression  to  a tropical  storm.  Debra  did  not  intensify  to  a 
tropical  storm  until  1800Z  on  28  August  1978.  The  WARF  position  esti- 
mates and  the  official  National  Hurricane  Center  (NHC)  track  are  shown 
in  Figure  2.  There  are  large  differences  between  the  WARF  position 
estimates  and  interpolated  positions  on  the  official  track  between 
1800Z  on  27  August  1978  and  0300Z  on  28  August  1978.  WARF  positions 
are  about  125  km  west  of  the  NHC  positions.  Part  of  this  difference 
can  be  explained  by  the  inability  by  all  sensors  to  precisely  locate 
the  storm  from  its  large  asymmetrical  center.  The  WARF  wind  map  shown 
in  Figure  la  is  representative  of  the  four  WARF  wind  maps  taken  at  this 
time. 

The  potential  sources  of  error  in  the  WARF-derived  positions  are 
associated  with  the  absolute  position  (range  and  azimuth)  accuracy  of 
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WARF-OERIVED  SURFACE  WIND  MAPS  FOR  DEBRA  MADE  AT  (a)  1801Z 
ON  AUGUST  27,  1978,  (b)  1600Z  ON  AUGUST  28,  1978,  (c)  2049Z  ON 
AUGUST  28,  1978,  AND  (d)  001 7Z  ON  AUGUST  29,  1978 


LATITUDE 


LONGITUDE  - 0 W 


FIGURE  1 (Continued) 


the  radar  and  the  location  of  a storm  center  from  the  wind  map.  The 
absolute  range  position  error  is  a function  of  the  midpath  ionospheric 
height.  Overhead  vertical-incidence  soundings  taken  at  WARF  generally 
provide  reasonable  estimates  of  the  ionospheric  height  at  midpath. 
Sometimes  large  position  errors  can  result  if  the  radar  signals  are 
actually  propagating  from  a different  ionospheric  layer  than  the  one 
determined  from  the  vertical- incidence  soundings.  We  minimized  this 
source  of  error  by  taking  radar  measurements  along  the  Mexican-United 
States  coastline.  The  land-echo  can  be  easily  distinguished  from  the 
sea-echo.  Since  the  coasL  is  approximately  perpendicular  to  the  radar 
beam,  the  location  of  the  land-sea  boundary  can  be  identified  in  range 
but  not  in  azimuth.  Using  this  method,  we  could  accurately  calculate 
the  ionospheric  height  at  midpath,  and  increase  the  accuracy  and 
reliability  of  the  range  position  estimates,  The  absolute  range  posi- 
tion estimate  is  no  more  than  10  km.  Azimuth  errors,  caused  by  iono- 
speric  tilting,  are  also  possible.  No  attempt  was  made  to  correct  for 
this  type  or  error.  Generally,  azimuth  errors  are  no  larger  than  25  km. 
The  surface  circulation  is  well  identified.  Positive  identification 
of  the  westerly  flow  north  of  the  storm  and  the  easterly  flow  south  of 
the  storm  indicate  that  the  storm  could  be  as  large  as  120  km  along  an 
axis  running  northeast.  We  placed  the  storm  center  along  an  axis  running 
parallel  and  halfway  to  these  winds.  It  is  usually  more  difficult  to 
identify  the  storm  center  in  range  because  of  the  left-right  ambiguity 
associated  with  each  wind  direction  measurement.  However,  the  northerly 
flow  emanating  from  the  NW  corner  of  the  map  and  the  southeasterly  flow 
emanating  from  the  SE  corner  of  the  map  are  unmistakable . The  size  of 
the  storm  center,  based  solely  on  the  WARF  wind  direction,  is  120  km 
x 120  km.  We  could  only  justify  the  selection  of  a storm  center  50  km 
further  out  in  range  based  on  surface  circulation  and  another  20  km 
based  on  the  absolute  position  error.  However,  we  would  still  be  more 
than  50  km  west  of  the  official  track  position.  The  consistency  of 
the  four  WARF  positions  and  the  high  quality  of  the  radar-wind  maps 
suggest  that  the  official  track  should  be  adjusted  westward. 


The  WARF-wind  maps  compiled  between  1000Z  on  28  August  1978  and 
0017Z  on  29  August  1978  are  shown  in  Figure  1(b)  through  (d).  The 
storm  intensified  during  this  period.  The  center  is  asymmetrical  with 
the  longest  axis  elongated  along  an  east-west  axis.  A reconnaissance 
aircraft  flight  through  the  center  at  1700Z  on  28  August  1978  measured 
the  center  as  75  km  x 150  km.  Peak  surface  winds  of  21  m/s  were  esti- 
mated from  the  aircraft. 

In  Table  1 we  compared  the  WARF  position  estimates  to  the  position 
estimates  made  from  the  reconnaissance  aircraft,  shore-based  Doppler 
microwave  radar,  and  satellite  cloud  photography.  The  position  esti- 
mates from  all  four  sensors  are  shown  in  Figure  3.  We  computed  the 
deviation  of  the  reconnaissance  aircraft  position  estimates  of  the 
storm  location  and  the  WARF,  satellite,  and  Galveston  radar  positions 
for  two  time  periods  at  approximately  2100Z  and  2400Z  on  28  August 
1978.  The  deviations  are  given  in  Table  2 for  2100Z  and  2400Z.  As 
expected,  the  agreement  between  the  Galveston  radar  and  the  aircraft 
was  the  best.  The  WARF  position  estimates  were  generally  south  of  the 
aircraft  fixes  and  the  satellite  positions  were  generally  north  of  the 
aircraft  fixes.  In  general  there  is  poor  agreement  between  all  four 
sensors  and  the  official  best  track  position  at  2100Z.  An  analysis 
and  discussion  of  the  deviation  of  the  satellite  and  aircraft  recon- 
naissance positions  from  the  official  track  is  given  in  Sheets  and 
Grieman.  The  average  deviation  of  the  satellite  positions  from  the 
official  track  for  data  analyzed  from  three  different  satellites  is 
about  65  km  but  this  average  deviation  decreases  for  intense  storms. 
The  average  deviation  of  the  aircraft  reconnaissance  fixes  from  the 
official  best  track  is  approximately  25  km.  Averaging  the  satellite 
and  skywave  radar  measurements  results  in  improved  position  accuracy 
and  takes  advantage  of  the  strengths  of  both  of  these  remote  sensing 
tools.  We  averaged  the  WARF  and  satellite  position  estimates  and  com- 
pared these  positions  to  the  aircraft  positions;  agreement  was  within 
20+5  km. 
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Table  1 

DEVIATION  OF  THE  WARF  RADAR  POSITION  ESTIMATE  FROM  THE 
RECONNAISSANCE  AIRCRAFT,  MICROWAVE  DOPPLER  RADAR,  AND 
SATELLITE  POSITION  ESTIMATES  FOR  TROPICAL  STORM  DEBRA 
AT  2100Z  AND  2400Z  AUGUST  28,  1978. 


Time 

Sensor 

Deviation 

(GMT) 

(km) 

2100 

Aircraft 

10 

Doppler  radar 

13 

Satellite 

53 

2400 

Aircraft 

40 

Doppler  radar 

32 

Satellite 

73 

Table  2 

DEVIATION  OF  THE  WARF  RADAR,  MICROWAVE  DOPPLER  RADAR, 
SATELLITE  CLOUD  PHOTOGRAPH,  WARF-SATELLITE  POSITION 
ESTIMATES  FROM  THE  RECONNAISSANCE  AIRCRAFT  POSITION 


i- 


ESTIMATES  FOR  TROPICAL 
AUGUST  28,  1978. 

STORM  DEBRA  AT  2100Z 

AND  2400Z 

Time 

(GMT) 

Sensor 

Deviation 

(km) 

2100 

WARF 

10 

. 

Satellite 

50 

Doppler  radar 

4 

WARF  + satellite 

25 

2400 

WARF 

40 

Satellite 

38 

Doppler  radar 

10 

WARF  + satellite 

15 
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FIGURE  3 COMPARISON  OF  THE  WARF-DERIVED  POSITIONS  WITH  THE  RECONNAISSANCE 
AIRCRAFT-,  SHORE-BASED  MICROWAVE  DOPPLER  RADAR-,  AND  SATELLITE- 
DERIVED  POSITIONS 
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IV  CONCLUSION 


We  have  compared  the  WARF  radar  measurements  to  the  official  NHC 
track  and  to  position  fixes  determined  by  other  conventional  sensors. 

We  attribute  the  large  difference  between  the  WARF  position  fixes  and 
the  NHC  track  before  1800Z  on  28  August  1978  to  an  eastward  bias  in  the 
NHC  track  caused  by  the  large,  poorly-defined  storm  center.  The 
agreement  between  the  WARF  position  fixes  and  the  reconnaissance  air- 
craft position  fixes  after  1800Z  was  within  25±15  km.  This  is  quite 
reasonable  considering  the  large  size  of  the  storm  center  associated 
with  tropical  storm  Debra. 

The  capability  to  produce  accurate  surface  wind  direction  maps  at 
the  time  and  location  specified  by  the  radar  operator  is  an  advantage  of 
skywave  radar  over  other  available  remote  sensing  techniques.  As  with 
early  satellite  cloud  photographs,  much  information  can  be  learned  as 
the  data  base  becomes  larger.  HF  skywave  radar  is  proving  itself  a 
valuable  and  reliable  remote  sensing  tool  to  monitor  tropical  storms. 
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HF  SKYWAVE  RADAR  MEASUREMENT 
OF  HURRICANE  WINDS  AND  WAVES 


Joseph  W.  Maresca,  Jr.  and  Christopher  T.  Carlson 


I.  INTRODUCTION 

We  measured  significant  wave  height,  and  surface  wind  speed  and 
direction  for  the  first  two  Gulf  of  Mexico  hurricanes  of  the  1977  season 
using  a high  frequency  (HF)  skywave  radar.  The  radar  measurements  were 
made  from  California  by  using  the  SRI-operated  Wide  Aperture  Research 
Facility  (WARF) . We  recorded  sea  backscatter  for  hurricanes  Anita  and 
Babe,  at  distances  more  than  3000  km  from  the  WARF,  by  means  of  single 
F-layer  ionospheric  reflection.  We  compiled  real-time  maps  of  the  sur- 
face wind  direction  field  within  a radial  distance  of  200  km  of  the 
storm  center,  then  estimated  the  hurricane  position  from  these  radar 
wind  maps,  and  developed  a track  for  Anita  over  a 4 -day  period  between 
30  August  and  2 September  1977  as  the  storm  moved  westward  across  the 
Gulf  of  Mexico.  The  radar  track  was  computed  from  17  independent  posi- 
tion estimates  made  before  Anita  crossed  the  Mexican  coast,  and  was 
subsequently  compared  to  the  official  track  produced  by  National  Hurri- 
cane Center  (NHC) . Agreement  between  the  WARF  position  estimates  and 
coincident  temporal  positions  on  the  NHC  smooth  track  was  ±19  km.  At 
aporoximately  0000Z  on  1 September  1977,  Anita  passed  within  50  km  of 
the  National  Data  Buoy  Off.'ce  (NDBO)  open  ocean  moored  buoy  EB-71,  and 
provided  us  with  the  opportunity  to  compare  WARF  estimates  of  the  sig- 
nificant wave  height,  and  surface  wind  speed  and  direction  in  all  four 
quadrants  of  the  storm  with  those  made  at  the  buoy.  Agreement  between 
the  WARF  and  EB-71  measurements  was  within  10%. 

Two  days  after  Anita  crossed  land,  tropical  storm  Babe — a weaker, 
short-lived  storm — developed.  WARF  estimates  of  the  significant  wave 
height,  and  surface  wind  speed  and  direction  were  made  for  selected 
regions  of  the  storm. 1 No  in  situ  wave  measurements  were  available  for 
comparison  to  the  WARF  measurements.  WARF  estimates  of  the  wind  speed 
were  compared  to  wind  speed  measurements  made  at  nearby  oil  platforms, 
and  surface  wind  speeds  computed  from  flight  level  winds  (305  m)  mea- 
sured by  a NOAA  reconnaissance  aircraft.  Agreement  was  again  within 
10%.  The  purpose  of  this  paper  is  to  describe  the  capability  of  re- 
motely monitoring  hurricanes  and  other  open  o_  an  storms  by  using  an  HF 
skywave  radar.  We  will  describe  the  important  aspects  of  the  WARF  sky- 
wave  radar,  the  sea-echo  Doppler  spectra,  the  method  of  analysis  used 
to  estimate  the  wave  and  wind  parameters,  and  the  accuracy  of  these 
radar-derived  quantities. 


SRI  International, 


Menlo  Park,  California  94025. 


II.  WARF  SKYWAVE  RADAR 


The  Wide  Aperture  Research  Facility  (WARF)2  is  a high-resolution, 
experimental, high-frequency  skywave  radar  located  in  central  California. 
The  radar  is  bistatic  and  operates  in  the  HF  band  between  6 and  30  MHz. 
Ocean  areas  are  illuminated  by  a 20-kW  swept-frequency  continuous-wave 
(SFCW)  signal  from  a transmitter  site  located  at  Lost  Hills,  California. 
The  energy  reflected  from  the  surface  beam  is  received  185  km  to  the 
north  at  Los  Banos,  California.  The  receiving  antenna  array  is  2.5-km 
long  and  consists  of  a double  linear  array  of  256  whip  antennas  produc- 
ing a nominal  1/2°  azimuthal  beamwldth  at  15  MHz.  The  signal  propagates 
to  and  from  remote  ocean  patches  by  means  of  one  or  more  ionospheric 
"reflections." 

The  WARF  coverage  area  is  shown  in  Figure  1.  The  radar  can  be  direct- 
ed either  east  or  west,  and  can  be  electronically  steered  in  azimuth  ±32° 
from  boresight  anywhere  within  the  coverage  area  in  1/4°  increments. 
Position  accuracy  is  a function  of  midpath  ionospheric  height  estimates 
whose  uncertainty  in  the  midpath  height  results  in  a nominal  position 
accuracy  of  approximately  20  km.  At  any  one  location,  the  accuracy  be- 
tween consecutive  measurements  in  range  and  azimuth  is  an  order  of  mag- 
nitude better.  WARF  has  multiple-beam  capability,  and  sea  backscatter 
is  usually  received  simultaneously  at  four  adjacent  ocean  areas  from 
four  different  beams  separated  by  1/4°.  The  size  of  the  ocean  scatter- 
ing patch  is  a function  of  the  beamwldth,  the  range,  the  range  cell 
separation,  and  the  number  of  range  cells  averaged  together.  The  size 
of  the  minimum  scattering  patch  at  a range  of  approximately  2000  km  is 
3 km  in  range  by  15  km  in  azimuth. 


III.  IONOSPHERIC  PROPAGATION 

The  ionosphere  consists  of  ions  produced  in  the  earth's  atmosphere, 
primarily  by  solar  radiation.  Radio-wave  propagation  by  means  of  iono- 
spheric reflection  occurs  primarily  between  elevations  of  100  km  and 
500  km.  A graph  of  electron  density  as  a function  of  height  may  show 
peaks  in  the  ionospheric  profile.  These  peaks  are  defined  as  layers 
and  are  designated  by  Eg  (sporadic-E) , E,  FI,  and  F2.  They  correspond 
to  peak  electron  densities  located  at  about  110,  120,  200,  and  300  km 
above  the  earth,  respectively.  Ionospheric  conditions  are  transient  in 
time  and  space,  and  they  depend  on  the  stability  and  strength  of  the 
electron  density  profile. 

The  minimum  radar  range  for  one  hop  ionospheric  propagation  is  ap- 
proximately 1000  km;  the  maximum  radar  range  is  approximately  3000  km. 
The  ionosphere  will  support  propagation  to  a specific  range  over  a 
limited  frequency  band.  The  achievable  range  depends  on  the  time  of 
day,  geographical  region,  and  ionospheric  height.  We  use  two  different 
types  of  real-time  ionospheric  soundings  at  WARF  to  manage  the  iono- 
spheric propagation.  An  oblique-incidence  sounding,  shown  in  Figure  2, 
is  primarily  used  to  determine:  the  relative  signal  strength;  the  radio 
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FIGURE  1 COVERAGE  AREA  OF  THE  WARF  HF  SKYWAVE  RADAR.  All  Anita 
measuicments  were  made  west  of  88°W  in  the  Gulf  of  Mexico. 


RADAR  FREQUENCY  - MHz 


frequencies  that  will  propagate  to  a particular  range;  and  certain  types 
of  Ionospheric  disturbances  such  as  traveling  waves,  and  focusing  or 
defocusing  of  energy.  A vertical-incidence  sound ingt shown  in  Figure  3, 
taken  between  the  WARF  transmitting  and  receiving  arrays  is  primarily 
used  to  measure  the  overhead  ionospheric  mode  structure  and  height  of 
each  ionospheric  layer.  A frequency  surveillance  spectrum  analyzer  is 
used  to  select  interference-free  frequency  bands. 


IV.  SEA  ECHO  DOPPLER  SPECTRUM 

The  sea  backscatter  received  at  the  WARF  is  coherently  processed  in 
range  and  Doppler  to  produce  a sea -echo  Doppler  spectrum.  We  usually 
process  21  independent  Doppler  spectra  spaced  at  3-km  range  intervals. 
These  spectra  are  obtained  simultaneously  at  each  of  four  adjacent  radar 
beams.  A total  of  84  independent  Doppler  spectra  are  obtained  for  each 
coherent  time  period.  We  compute  an  average  spectrum  from  a subset  of 
these  Doppler  spectra,  depending  on  the  type  of  measurement  and  the  time 
and  space  scales  associated  with  the  ocean  surface  features.  An  example 
of  a mean  sea -echo  Doppler  spectrum  produced  by  averaging  112  spectra 
obtained  from  four  consecutive  102.4  s coherent  time  periods,  over  a 
scattering  patch  consisting  of  21  range  cells  and  3 adjacent  beams, is 
shown  in  Figure  4. 

The  sea-echo  Doppler  spectrum  shown  in  Figure  4 is  characterized 
by  two  dominant  first-order  echoes  surrounded  by  a second-order  con- 
tinuum. Crombie**  Interpreted  the  first-order  echoes  in  terms  of  simple 
Bragg  scattering  that  represented  a resonant  response  between  radio 
waves  of  wave  number  kQ  and  ocean  waves  of  wave  number  k = 2 k0.  The 
radar  measures  the  relative  power  and  Doppler  of  the  ocean  waves  travel- 
ing radially  toward  or  away  from  the  radar.  The  power  ratio  of  the  two 
first-order  echoes  are  indicative  of  the  wave  direction  of  the  waves  of 
wave  number  k.  Because  k is  usually  large  (k  > 0.5),  it  is  assumed  that 
the  wind  direction  is  identical  to  the  direction  of  these  waves. 

The  wave -height  spectrum  is  derived  from  the  second-order  structure 
surrounding  the  first-order  echoes.  For  hurricanes,  the  power  in  the 
second-order  echoes  is  large.  As  the  total  wave  energy  increases,  the 
anplitude  of  the  second-order  echoes  increases  as  illustrated  in 
Figure  5.  Barrlck^*^  derived  theoretical  expressions  that  accurately 
model  the  HF  scattering  process  to  second  order.  For  a specific  direc- 
tional wave  spectrum,  the  model  computes  the  Doppler  spectrum.  The  ef- 
fects of  the  wind  direction,  wave  directionality  and  the  wave  frequency 
spectrum  on  the  modeled  Doppler  spectrum  have  been  extensively  studied 
through  the  use  of  this  model. 


V.  HURRICANE  DATA  SAMPLING 

Data  sampling  during  a hurricane  is  divided  into  two  tasks  to  opti- 
mize the  sampling  time  and  the  data  quality.  The  spectral  resolution. 


VERTICAL  FREQUENCY  - MHz 

FIGURE  3 VERTICAL-INCIDENCE  IONOSPHERIC  SOUNDING  TAKEN  AT  WARF. 

The  vertical-incidence  sounding  is  a plot  of  overhead  energy  for  different 
frequencies. 


FIGURE  4 AVERAGE  SEA  ECHO  DOPPLER  SPECTRUM  RECORDED  WITHIN  35  km 

OF  THE  CENTER  OF  HURRICANE  ANITA  AT  2343Z  ON  AUGUST  31,  1977 
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FIGURE  5 EXAMPLE  OF  TWO  SYNTHETIC  DOPPLER  SPECTRA  (b)  PRODUCED  FROM 


TWO  INPUT  WAVE  SPECTRA  (a)  WITH  THE  SAME  DIRECTIONAL  DISTRIBUTION 
AND  RADAR-TO-WIND  DIRECTION,  BUT  DIFFERENT  TOTAL  WAVE  ENERGY 
(0.02  Hz  resolution) 


directly  related  to  the  coherent  integration  time,  can  be  much  coarser 
for  first-order  measurements  than  for  second-order  measurements.  Wind 
direction  estimates  are  computed  from  the  first-order  echoes,  and  can 
be  computed  considerably  more  rapidly  than  wave  height  and  wind  speed 
estimates,  which  are  computed  from  the  second-order  echoes.  Usually, 
the  longer  the  coherent  integration  time,  the  greater  the  influence 
the  ionosphere  has  on  the  quality  of  the  data. 

The  quality  of  the  recorded  sea  backscatter  depends  on  the  iono- 
spheric conditions  over  short  periods — on  the  order  of  minutes.  High- 
quality  sea  backscatter  is  obtained  if  the  radio  waves  propagate  by 
means  of  a strong,  single,  stable,  coherent  ionospheric  layer.  Some- 
times the  signals  may  be  received  at  the  same  time  from  two  or  more 
different  paths  (multipath).  In  this  case,  the  second  or  succeeding 
signals  will  be  reflected  from  different  parts  of  the  ocean  and  differ- 
ent parts  of  the  ionosphere,  and  will  contaminate  the  sea  echo  received 
from  the  first  path.  If  the  ionosphere  is  changing  in  time  or  space 
during  the  coherent  radar  dwell  (time  period),  further  degradation  of 
the  data  will  occur.  The  ability  to  predict  the  ionospheric  conditions 
would  enable  the  radar  operator  to  minimize  the  contaminating  effects 
of  the  ionosphere,  improve  the  quality  of  sea  backscatter,  and  reduce 
the  sampling  time.  The  ionospheric  soundings  provide  some  data  quality 
Information.  The  vertical  and  oblique  incidence  soundings  are  taken 
every  10  minutes;  a complete  sounding  requires  approximately  3 minutes. 
The  coherent  radar  measurements  made  at  WARF  require  between  10  and 
100  seconds  to  complete.  Because  the  time  required  to  complete  a sound- 
ing is  greater  than  the  time  required  to  record  the  sea  backscatter 
data,  assessment  of  the  data  quality  is  difficult  for  rapidly  changing 
ionospheres.  Therefore,  real-time  output  of  the  data  from  the  WARF 
site  minicomputer  is  used  to  verify  data  quality. 

The  wind  direction  measurement  is  not  extremely  sensitive  to  iono- 
spheric contamination  caused  by  multipath  or  smearing  because  only  the 
amplitude  of  the  two  etrong  first-order  echoes  must  be  measured.  A coher- 
ent integration  time  of  12.8  seconds  (0.078-Hz  resolution)  is  suffi- 
cient to  resolve  the  peaks  of  the  first-order  echoes.  We  can  map  the 
wind-direction  field  in  a hurricane  by  scanning  in  range  and  azimuth. 

It  is  possible  to  routinely  map  the  surface-vlnd-direction  field  of  a 
hurricane  and  this  can  be  accomplished  in  about  10  minutes.  Once  the 
surface-wind-direction  map  Is  made,  the  storm  center  can  be  Identified 
for  tracking  purposes,  and  regions  of  interest  can  be  selected  for  more 
extensive  monitoring  of  wind  speed  and  wave  height  anywhere  within  the 
storm. 

The  significant  wave  height  and  wind  speed  measurements  are  sensi- 
tive to  ionospheric  contamination.  This  contamination  is  the  largest 
source  of  error  in  these  measurements.  A coherent  integration  time  of 
102.4  seconds  (0.01-Hz  resolution)  is  required  to  resolve  the  second- 
order  echoes.  The  ionosphere  does  not  generally  support  coherent  inte- 
gration time  periods  of  this  length.  Multipath  and  ionospheric  smearing 
can  seriously  degrade  the  weaker  second-order  echoes.  Because  of  this 
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contamination,  we  are  not  able  to  routinely  estimate  wave  height  for 
each  102.4-second  time  period  as  we  are  able  to  calculate  wind  direc- 
tion for  each  12.8-second  time  period.  A sampling  strategy  that  com- 
bines careful  propagation  management  through  selection  of  frequencies, 
which  result  in  a stable,  coherent,  single  propagation  path,  and  signal 
processing  that  minimizes  the  contaminating  effects  of  the  ionosphere 
are  used  to  obtain  a data  set  suitable  for  analysis.  Recent  work  by 
SRI  and  NOAA^  has  resulted  in  improved  methods  of  collecting  high-qual- 
ity data  by  sorting  the  data  according  to  a spectral  sharpness  index. 

The  effect  of  ionospheric  contamination,  however,  is  less  severe  for  data 
recorded  during  large  waves  generated  during  a hurricane.  The  amplitude 
of  the  second-order  echoes  containing  the  wave  height  Information  may  be 
stronger  than  the  contamination  effects,  and  thus,  wave  height  can  be 
calculated  despite  the  contamination. 

For  the  Anita  and  Babe  wind  direction  measurements,  we  divided  the 
data  into  16  groups  and  analyzed  three  consecutive  12.8-second  coherent 
radar  dwells.  Each  wind  direction  estimate  was  calculated  from  a mini- 
mum of  15  Doppler  spectra.  At  a range  of  3000  km,  the  size  of  each 
scattering  patch  was  15  km  x 25  km.  It  would  be  desirable  to  compute 
wave  height  and  wind  speed  from  a similar  data  set,  but  this  is  not 
generally  possible.  Longer  coherent  time  periods  and  more  independent 
samples  of  the  spectra  are  required  to  obtain  a high  quality  sample. 

We  could  collect  the  data  over  a small  scattering  patch  by  averaging 
over  a long  time,  or  we  could  increase  the  scattering  patch  size  and 
average  in  space.  Averaging  in  space  is  preferable  because  it  reduces 
the  total  time  required  to  obtain  a mean  Doppler  spectrum.  For  the 
Anita  and  Babe  wave  height  and  wind  speed  measurements,  we  analyzed  the 
data  from  three  of  the  adjacent  azimuth  cells  and  21  contiguous  range 
cells.  The  total  scattering  patch  was  63  km  x 50  km.  Several  consecu- 
tive integration  periods  are  required  to  record  the  data. 


VI.  WIND  DIRECTION  MEASUREMENT 

HF  skywave  radar  has  been  used  to  map  the  surface-wind  fields  asso- 
ciated with  large  weather  systems®  and  tropical  storms. 9 The  radar 
measured  surface-wind  directions  are  derived  from  the  predominant  direc- 
tion of  ocean  gravity  waves,  approximately  10-m  long.  The  waves  satis- 
fying the  first-order  Bragg  scattering  condition,  k - 2k0,  are  assumed 
to  be  tightly  coupled  to  the  wind  for  time  scales  on  the  order  of  tens 
of  minutes.  This  assumption  is  reasonable  for  the  high  wind  speed  con- 
ditions associated  with  hurricanes.  Available  directional  wave  spectra 
measurements 10~12  indicate  that  the  dominant  wave  direction  is  repre- 
sentative of  the  predominant  wind  direction.  For  open  ocean  conditions, 
agreement  between  the  WARF  radar  and  shipboard  anemometer  measurements 
of  wind  direction  is  ±16°.  For  hurricane  winds,  the  agreement  be- 
tween coincident  wind  direction  measurements  made  by  the  NOAA  National 
Data  Buoy  Office  (NDBO)  data  buoys  and  the  WARF  radar  is  better  than 
10°.l»* 
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The  radar  measures  the  relative  power  between  the  approaching  and 
receding  waves  that  satisfies  the  Bragg  scattering  condition.  If  a 
cosine  directional  distribution*** 

G(0)  - coss0/2  (1) 

is  assumed,  then  the  relative  power  of  the  approaching  and  receding 
waves  measured  by  the  radar  is  sufficient  to  estimate  0 with  an  ambi- 
guity about  the  beam  direction.  This  left-right  ambiguity  is  resolved 
by  the  predictable  cyclonic  surface  circulation  within  the  hurricane. 

The  shape  of  G(0)  is  controlled  by  the  spreading  parameter,  s,  where  0 
is  the  angle  between  the  radar  beam  and  the  wind  direction.  For  open 
ocean  conditions,  we  have  estimated  s from  several  models. **» 13  For 
the  maximum  hurricane  winds,  the  values  of  s estimated  from  these  models 
are  too  low.  Based  upon  previous  hurricane  analyses  and  spot  measure- 
ments of  wind  direction  at  NDBO  data  buoys,  we  used  values  of  s between 
1.0  and  2.0.  No  attempts  were  made  to  account  for  variations  in  s as  a 
function  of  location  within  the  hurricane. 


VII.  SIGNIFICANT  WAVE  HEIGHT 

Barrick^  derived  an  integral  expression  that  predicts  the  Doppler 
spectrum  for  a specific  directional  wave  spectrum  input.  Recent  ef- 
forts have  succeeded  in  inverting  this  Integral  expression  to  compute 
the  input  rms  wave  height, *^»*^  one-dimensional  wave  frequency  spec- 
trum, lo-19  aiuj  t^e  directional  distribution.^-^  Barrick' s*^» *^  ex- 
pressions' have  been  used  to  analyze  skyvave  radar  data  recorded  for  a 
Pacific  Ocean  storm*-*  and  tropical  storms.  1*20 

We  used  a power  law  derived  from  simulated  data  by  Maresca  and 
Georges *5  to  compute  rms  wave  height  by  relating  the  ratio  of  the  total 
second-order  and  first-order  power  to  the  rms  wave  height: 

kQh  - aR*  (2) 

where  0.2  <-k-h  < 1.0,  h is  the  rms  wave  height;  Icq  is  the  radar  wave 
number;  R2  is  the  ratio  of  the  total  second-  to  total  first-order  power; 
and  a ■ 0.8  and  b * 0.6  are  constants.  This  average  expression  was  de- 
rived from  theoretical  simulations  of  the  Doppler  spectra  for  different 
radar-to-wind  directions,  directional  distributions,  functional  forms 
of  the  wave- frequency  spectrum,  and  operating  radar  frequencies.  Equa- 
tion (2)  is  accurate  to  within  10%.  Discussion  of  the  errors  can  be 
found  in  Maresca  and  Carlson,*  and  Maresca  and  Georges.*** 


VIII.  WIND  SPEED 


Historically,  wave  models  have  been  developed  to  predict  wave 
height  and  the  wave  spectrum  from  an  Input  wind  field.  The  accuracy 
of  these  models  Is  dependent  upon  the  accuracy  of  the  Input  winds. 
Hasselmann  et  al.^l  proposed  a one-dimensional  parametric  wind-wave 
model  for  fetch  limited  growing  wind-sea  conditions.  Ross  and 
Cardone22-25  empirically  derived  a power-law  expression  for  hurricanes 
based  on  the  form  proposed  by  Hasselmann  that  relates  the  nondlmensional 
wave  energy,  E,  by  using  wind,  wave, and  fetch  measured  during  hurricanes 
Ava,  Camille,  and  Eloise.  For  hurricanes, 

E - 2.5  x 10-5  R0,45  (3) 


where  E * Eg^/W*;  R - rg/W^;  E ■ h^;  and  Hs  ■ 4h.  In  E and  R,  E is  the 
total  wave  energy;  h is  the  rms  wave  height;  Hs  is  the  significant  wave 
height;  r is  the  radial  distance  from  the  eye.  to  the  measurement  point 
that  accounts  for  fetch;  g is  the  gravitational  acceleration;  and  W is 
the  wind  speed.  Solving  for  wind  speed  in  Eq.  (3),  we  obtain 
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(4) 


The  wind-wave  model  used  to  derive  Eq.  (4)  is  applicable  for  slow  moving 
storms  in  which  W > 15  m/s  and  R < 3 x 10*.  For  the  unusual  cases  where 
the  storms  move  very  fast  or  very  slow,  Ross  and  Cardone^*  showed  that 
significant  differences  in  the  modeled  and  measured  wave  heights  occur. 


We  used  Eq.  (4)  to  calculate  wind  speed  for  both  Anita  and  Babe 
and  compared  our  results  with  the  wind  speeds  measured  at  NDBO  buoys, 
oil  platforms,  and  by  reconnaissance  aircraft.  The  radial  fetch  (r) 
was  measured  from  the  WARF-derived  wind  maps,  and  the  wave  height  (h) 
was  computed  using  Eq.  (2).  The  radar-derived  W is  not  an  instantane- 
ous wind  speed  estimate;  it  is  a smooth  temporal  and  spatial  average 
of  the  winds.  Our  radar-derived  W was  compared  to  the  15-minute  wind 
speed  averages  made  at  NDBO  moored  data  buoys. 


IX.  ANITA  MEASUREMENTS 

Hurricane  Anita  formed  as  a tropical  depression  in  the  Gulf  of 
Mexico  at  about  1200Z  on  29  August  1977 . Anita  developed  into  a tropi- 
cal storm  at  approximately  0600Z  on  30  August  1977,  and  about  12  hours 
later  intensified  into  the  first  Gulf  of  Mexico  hurricane  of  the  1977 
season.  As  Anita  moved  west  across  the  Gulf,  winds  in  excess  of  75  m/s 
were  recorded.  Five  days  of  skywave  data  beginning  29  August  1977, 
were  recorded  prior  to  Anita's  landfall,  2 September  1977,  approxi- 
mately 48  km  south  of  Brownsville,  Texas.  Twenty-one  radar  wind  maps 
were  compiled  at  WARF.  The  first  4 wind  maps  were  not  used  in  the 
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radar-derived  track  presented  here  because  the  radar  showed  two  distinct 
centers  during  this  early  period.  On  30  August  1977,  the  storm  intensi- 
fied and  developed  one  center.  The  wind  maps  were  updated  3 to  5 times 
per  day  during  both  daytime  and  nighttime  periods  and  were  used  to  de- 
velop the  WARF-derived  track.  Figure  6 shows  the  radar-derived  posi- 
tions in  relation  to  the  official  NHC  smooth  track  produced  from  recon- 
naissance aircraft  measurements,  visible  and  infrared  satellite  cloud 
photographs,  and  shore-based  microwave  Doppler  radar.  The  relative 
agreement  between  the  WARF  position  estimates  and  the  interpolated 
temporal  position  estimates  along  the  smooth  track  is  ±19  km. 

There  are  two  potential  sources  of  error  associated  with  the  WARF 
hurricane  position  fixes:  the  absolute  position  error  of  the  radar  con- 
sisting of  range  and  azimuth  errors,  and  the  errors  associated  with  lo- 
cating the  storm  center  from  the  radar  wind  direction  measurements.  We 
estimate  the  range  errors  of  the  radar  caused  by  errors  in  determining 
the  ionospheric  height  at  midpath  to  be  20  km.  If  a coastal  scan  is  in- 
cluded as  part  of  collecting  the  wind  map  data,  the  land  echo  can  be  used 
.as  a reference  to  more  accurately  determine  the  ionospheric  height,  and 
therefore,  reduce  this  error.  We  estimate  the  error  in  azimuth  caused 
by  ionospheric  tilting  to  be  20  km.  These  range  and  azimuth  errors  can 
be  reduced  significantly  by  installing  an  HF  repeater  along  the  coast  which 
receives  signals  and  transmits  them  back  with  a known  frequency  shift. 

When  we  assume  similar  mean  ionospheric  conditions  within  200  km  of  the 
storm  center,  the  entire  wind  map  can  be  translated  in  azimuth  and  range 
to  correct  for  the  absolute  position  error.  The  location  of  the  wind 
direction  measurement  with  respect  to  the  storm  center  is  generally  not 
affected  by  these  position  errors.  The  error  associated  with  determining 
the  storm  center  from  the  radar  maps  is  about  20  km.  The  error  is 
caused  by  the  left/right  ambiguity  In  the  wind  direction  measurement. 

The  average  maximum  error  from  these  two  potential  sources  of  error  is 
about  40  km.  In  comparing  the  WARF  position  fixes  to  the  NHC  track  we 
found  relative  differences  of  between  5 and  50  km,  and  these  relative 
differences  can  be  attributed  to  the  sources  of  error  just  discussed. 

Anita  passed  50  km  south  of  NDBO  buoy  EB-71  at  about  0000Z  on 
1 September  1977.  Two  WARF-derived  wind  maps  were  made  at  2140Z  on 
30  August  1977  and  0120Z  on  1 September  1977,  which  brackets  this  time 
period.  One  of  these  wind  maps  is  shown  in  Figure  7.  Also  shown  on 
Figure  7 are  surface  wind  direction  fields  derived  from  data  recorded 
by  NDBO  buoy  EB-71.  These  buoy-measured  wind  directions  were  recorded 
at  2-hour  intervals  during  the  period  ±18  hours  of  Anita's  passing 
EB-71.  The  buoy-derived  wind  field  was  computed  by  a time-space  conver- 
sion that  assumed  uniform  wind  direction  and  lateral  storm  motion  during 
this  period.  We  compared  the  buoy-derived  wind  directions  to  the  WARF- 
derived  wind  directions;  agreement  was  within  19°.  Agreement  between 
the  WARF-derived  wind  direction  estimate  coincident  in  rime  and 
space  with  the  buoy  wind  direction  estimate  was  1°. 

Between  2314Z  on  31  August  1977  and  0020Z  on  1 September  1977,  WARF 
measurements  were  made  at  five  locations  surrounding  the  center  of  the 
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FIGURE  7 WAHF-DERIVED  WIND  DIRECTION  MAP  MADE  FOR  ANITA  AT  2140Z  ON 
AUGUST  31,  1977 


storm.  The  location  of  each  measurement  relative  to  the  storm  center 
was  Interpolated  from  the  two  wind  maps.  We  computed  the  wind  direction 
by  using  Eq.  (1),  wind  speed  by  using  Eq.  (4),  and  wave  height  by  using 
Eq.  (2)  at  each  location  (see  Table  1),  and  compared  these  measurements 
to  a buoy-derived  wind  and  wave  field.  The  maps  of  the  spatial  distri- 
bution of  the  wind  direction,  wind  speed,  and  wave  height  were  compiled 
from  NDBO  EB-71  data  buoy  measurements.  Each  parameter  was  plotted  in 
relation  to  the  storm  center;  they  are  shown  in  Figures  8,  9,  and  10. 

We  assumed  that  Anita  moved  uniformly  with  no  change  in  the  meteorological 
conditions  during  the  period  18  hours  before  and  18  hours  after  passing 
the  buoy. 


Table  1 

WARF  ESTIMATES  OF  SIGNIFICANT  WAVE  HEIGHT, 
WIND  SPEED,  AND  WIND  DIRECTION 


Point 

Lat 

(°N) 

Long 

(°W) 

Time 

(GMT) 

' r*  ' 
(km) 

N* 

HS* 

(m) 

w*-" 

(m/s) 

(°N) 

A 

25.7 

92.9 

2314 

35 

80 

5.8 

26.7 

277.5 

A 

25.7 

92.9 

2343 

35 

112 

5.2 

22.8 

— 

B 

26.3 

92.1 

2324 

75 

80 

6.0 

24.4 

95.1 

C 

26.3 

93.1 

2358 

65 

35 

5.8 

24.4 

70.2 

D 

25.7 

92.1 

0003 

65 

134 

5.1 

22.5 

168.8 

E 

25.2 

91.1 

0020 

180 

49 

4.6 

18.1 

137.2 

r * Radial  Distance 
N = Number  of  Spectra  Averaged 
Hs=  Significant  Wave  Height 
W = Wind  Speed 
4>  ■ Wind  Direction 


Tne  significant  wave  height  shown  in  Figure  8 was  measured  at  the 
buoy  every  3 hours;  the  wind  direction  and  wind  speed  shown  in  Figures  9 
and  10  were  measured  at  the  buoy  every  2 hours.  During  this  36-hour 
period,  Anita  began  to  intensify,  and  the  validity  of  the  buoy-derived 
wind  and  wave  fields  are  suspect.  Exact  comparison  of  the  EB-71  and 
WARF  measurements  are  difficult  because  of  the  differences  in  the  time, 
location,  aH  area  of  ocean  monitored.  On  Figure  8 we  also  included 
the  wave  forecast  for  significant  wave  height  computed  by  Cardone  et 
al.  5 for  comparison. 

The  WARF  estimates  made  at  Point  B (26.3°N,  92.1°W)  were  in  close 
proaiaity  to  the  buoy-derived  estimates  located  at  (26.2°N,  92.1° W). 

The  remaining  WARF  wind  and  wave  height  estimates  were  too  far  away  from 
' H*  buoy-derived  quantities  for  direct  comparison,  but  the  agreement 
•M>'veen  the  WARF-  and  buoy-derived  wind  and  wave  fields  was  reasonable. 


88 


© BUOY-DERIVED  WAVE  HEIGHTS 
• WARF-DERIVED  WAVE  HEIGHTS 


FIGURE  8 COMPARISON  OF  THE  WARF-DERIVED  SIGNIFICANT  WAVE  HEIGHTS  (m) 
MEASURED  BETWEEN  2314Z  ON  AUGUST  31,  1977  AND  0020Z  ON 
SEPTEMBER  1,  1977  AND  THE  EB71-DERIVED  SIGNIFICANT  WAVE 
HEIGHTS  MADE  BETWEEN  0600Z  ON  AUGUST  31,  1977  AND  1800Z  ON 
SEPTEMBER  1,  1977.  The  wave  height  contours  are  reproduced  from 
Figure  9 of  Reference  25.  The  letter  designations  are  given  on  Figure  9. 
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© BUOY-DERIVED  WIND  DIRECTIONS 

• WARF-DERIVED  WIND  DIRECTIONS 
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LONGITUDE  - °W 


FIGURE  9 COMPARISON  OF  THE  WARF-DERIVED  WIND  DIRECTIONS  (-♦-)  MADE 

BETWEEN  2314Z  ON  AUGUST  31,  1977  AND  0020Z  ON  SEPTEMBER  1,  1977 
AND  THE  EB71-DERIVED  WIND  DIRECTIONS  (— — ) MADE  BETWEEN  0600Z 
ON  AUGUST  31,  1977  AND  1800Z  ON  SEPTEMBER  1,  1977 
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The  wave  forecast  was  compared  to  both  the  buoy-  and  WARF-derived  wave 
heights.  We  found  good  agreement  between  the  forecast,  buoy,  and  WARF 
estimates  of  wave  height  along  the  4.6-m  contour  east  of  the  hurricane. 
The  buoy  estimate  west  of  the  hurricane  along  the  4.6-m  contour  line 
suggests  the  radius  of  this  contour  line  was  too  large.  WARF  estimates 
of  wave  height  at  Points  B and  C,  located  at  the  extremes  of  the  fore- 
cast region  of  highest  waves,  were  lower  than  the  hindcast.  The  good 
agreement  of  Point  B with  the  nearby  buoy  estimate,  and  the  consistency 
of  the  WARF  estimates  at  Point  B and  C suggest  that  the  7.6-m  contour 
should  have  been  smaller.  The  WARF  estimates  of  wave  height  at  Points  A 
and  D are  located  between  the  6.1-m  contours.  Again,  the  buoy-  and  WARF- 
derived  wave  heights  suggested  that  the  6.1-m  contour  line  was  too  large. 
We  should  also  note  that  any  contamination  of  the  sea-echo  Doppler 
spectra  by  the  ionosphere  would  result  in  radar  wave-height  estimates 
that  would  be  too  high.  The  composite  of  wave  height  data  obtained 
from  the  hindcast  and  EB-71  buoy  indicated  the  validity  of  the  WARF 
wave  height  estimates. 

The  agreement  between  the  WARF-derived  estimates  of  wind  speed  and 
the  buoy-derived  estimates  of  wind  speed  is  good.  There  are  three 
principal  errors  associated  with  the  WARF  wind  speed  estimate:  error 
in  estimating  the  radial  fetch,  error  in  estimating  the  rms  wave  height, 
and  error  in  the  parametric  model.  We  computed  the  error  in  calculat- 
ing wind  speed  for  a ±0.5-m  error  in  estimating  wave  height  for  a signi- 
ficant wave  height  of  5.5  m (9.1%  error)  and  for  radial  fetches  of  30, 

50,  70,  and  100  km.  The  errors  were  less  than  1.6  m/s.  We  also  computed 
the  error  in  calculating  wind  speed  for  a ±20-km  error  in  estimating 
the  radial  fetch  for  a significant  wave  height  of  5.5  m and  radial 
fetches  of  30,  50,  70,  and  100  km.  For  radial  fetches  greater  than 
30  km,  a +20-km  error  causes  an  error  of  less  than  2 m/s  in  wind  speed. 
For  radial  fetches  greater  than  50  km,  a -20-km  error  causes  an  error 
of  less  than  2 m/s  in  wind  speed.  This  represents  less  than  an  8%  error. 
These  errors  are  typical  of  the  WARF  estimates  of  the  significant  wave 
height  and  radial  fetch  measurements.  The  errors  associated  with  the 
model  are  discussed  by  Ross  and  Cardone.24  For  Anita,  the  mean  and  rms 
differences  between  the  Cardone  et  al.^5  parametric  model  forecast  and 
measured  wave  heights  at  EB-71  is  0.21  ±0.83  m.  This  includes  errors 
in  measuring  wave  height  at  the  buoy  and  in  radial  fetch  from  the  con- 
ventional position  fixes.  We  also  calculated  the  wind  speed  using  Eq. 

(4)  for  some  of  the  buoy-measured  wave  heights  shown  in  Figure  8,  and 
compared  the  calculated  wind  speed  measurements  to  wind  speeds  measured 
at  the  buoy  (Table  2) . The  data  is  indicative  of  the  accuracy  we  could 
expect  from  the  WARF  estimates  of  wind  speed  using  Eq.  (4).  For  these 
data,  we  believe  the  largest  sources  of  error  in  the  comparison  were 
the  uncertainty  in  the  radial  distance  to  each  point  caused  by  compiling 
the  map  over  a 36-hour  period,  and  the  assumption  of  a symmetrical  dis- 
tribution of  the  winds. 
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Table  2 


COMPARISON  OF  WIND  SPEED  VALUES  CALCULATED 
FROM  EQ.  (4)  DERIVED  FROM  EB-71  SIGNIFICANT 
WAVE  HEIGHT  MEASUREMENTS 


Lat 

(°N) 

Long 

<°W) 

V 

(m) 

r* 

(km) 

u * 
(m?s) 

TJ  * 

wm 

(m/s) 

Difference 

(m/s) 

25.5 

94.8 

2.5 

217 

11.9 

7.6 

44.3 

25.7 

94.4 

2.9 

174 

13.5 

9.0 

44.5 

25.7 

93.9 

3.1 

124 

14.9 

13.3 

4-1.6 

25.8 

93.4 

5.5 

69 

23.4 

17.4 

+6.0 

26.0 

92.5 

6.5 

24 

30.4 

34.1 

-3.7 

26.3 

91.5 

4.7 

126 

19.4 

23.1 

-3.7 

26.6 

90.7 

4.6 

212 

17.7 

18.1 

-0.4 

Hs  « Significant  Wave  Height 
r = Radial  Distance 
Wc  = Computed  Wind  Speed 
Wm  = Measured  Wind  Speed 


X.  BABE  DATA 

We  recorded  skywave  data  for  Babe  from  2000Z  on  4 September  1977 
until  its  landfall  on  5 September  1977.  Wind  direction  maps  were  made  at 
2253Z  on  4 September  1977  and  1507Z  on  5 September  1977.  Babe  was  al- 
ready onshore  before  the  second  map  was  made.  We  analyzed  a sea  echo 
Doppler  spectrum  near  the  peak  winds  for  wind  speed.  This  spectrum 
was  recorded  at  a radius  of  50  km  from  the  hurricane  center  (27.9°N, 
91.6°W).  The  WARF  estimate  of  significant  wave  height  was  3.6  m.  The 
wind  speed  was  computed  by  using  Eq.  (4),  and  was  compared  to  wind  speeds 
measured  at  several  offshore  oil  platforms  at  2100Z  and  from  aircraft 
reconnaissance  at  1800Z.  Winds  measured  from  aircraft  at  an  altitude 
of  305  m were  reduced  to  the  equivalent  10  m wind  for  comparison  to  the 
radar  data.  We  reduced  the  aircraft  winds  to  the  surface  using  a simple 
ratio  relating  upper  level  gradient  wind  to  the  surface  wind,  indicated 
by  Elsberry  et  al.26  calculations.  Using  the  two  layer  Cardone27  marine 
boundary  layer  model,  Elsberry  et  al.26  computed  the  ratio  of  the  wind 
at  the  top  of  the  upper  layer  to  the  wind  at  the  top  of  the  surface 
layer  for  different  regions  of  the  hurricane,  different  surface  rough- 
ness and  different  ratios  of  heat  conductivity  to  eddy  viscosity.  For 
moderate  to  high  wind  speeds,  the  top  of  the  surface  layer  is  approxi- 
mately 20  m.  This  wind  ratio  ranges  from  about  0.5  to  0.85.  The  lower 
value  represents  regions  near  the  peak  winds.  We  assumed  the  305-m  air- 
craft wind  was  representative  of  the  wind  at  the  top  of  the  upper  layer 
and  reduced  it  to  20  m using  a ratio  of  0.6.  The  corrected  20-m  wind 
is  18.5  m/s.  We  realize  the  error  associated  with  this  calculation  can 


93 


be  large.  We  computed  the  19.5-m  WARF  wind  speed  for  comparison  by 
assuming  a logarithmic  profile  in  the  surface  boundary  layer  of  the 
form 


where  W(z)  is  the  wind  speed  at  a height  of  z,  W*  = CqW^  is  the  fric- 
tion velocity,  CD  is  the  drag  coefficient,  k 5 0.4  is  Von  Karman's  con- 
stant, and  zQ  is  the  surface  roughness.  We  estimated  z0  = 0.00392  m for 
hurricane  conditions  using  the  wind  speeds  measured  during  Eloise^®  at 
NDBO  buoy  EB-10  at  10  m,  and  a constant  drag  coefficient  of  0.0026  pro- 
posed by  Wu^9  for  high  wind  speeds  to  calculate  W*.  We  found  good 
agreement  between  our  19.5~m  computed  wind  speed  using  Eq.  (5)  for 
Eloise  and  the  19.5-m  Eloise  wind  speed  computed  by  Ross  and  Cardone.^ 
Using  zQ  = 0.00392  m and  Cp  = 0.0026  m,we  computed  the  19.5-m  WARF  Babe 
wind  speed  to  be  18.3  m/s.  We  also  computed  the  maximum  wind  speed, 
wmax*  anc*  the  maximum  sustained  wind  speed,  Ws,  for  a storm  moving  at 
5.1  m/s  from  the  following  expressions^  relating  the  central  pressure 
and  radius  of  the  storm  to  wind  speed  at  10  m: 

W = 0.868  [6.45 (P  -P  )1/2  - 0.296  rf]  (6) 

max  1 n o' 


and 


W = 0.865  W + 0.5  V_,  (7) 

s max  F 

where  ^max  and  Ws  are  in  m/s,  Pn  is  the  normal  pressure  of  1013  mbars,  PQ 
is  the  central  pressure  in  mb,  r is  the  radius  in  km,  f is  the  coriolis 
parameter  in  radians/hour,  and  Vp  is  the  forward  motion  of  the  storm  in 
m/s.  The  maximum  wind  recorded  from  the  aircraft  at  1800Z  were  located 
approximately  60  km  from  the  center  and  the  central  pressure  was  re- 
ported as  1000  mb.  We  calculated  * 16.4  m/s  and  Ws  = 16.8  m/s. 

We  summarize  these  wind  speed  estimates  in  Table  3.  Because  none  of 
the  wind  speed  estimates  were  coincident  in  time  or  space  with  the  WARF 
estimate  and  the  assumptions  inherent  in  deriving  these  quantities  we 
only  compared  the  results  qualitatively.  The  WARF  wind  speed  estimate 
is  quite  reasonable. 


XI.  SUMMARY 

Spatially-averaged  hurricane  wind  speed,  wind  direction,  and  wave 
height  estimates  made  at  the  WARF  for  Anita  and  Babe  were  compared  to 
point  measurements  made  at  NDBO  buoys  and  oil  platforms  and  by  recon- 
naissance aircraft.  Agreement  was  within  the  nominal  measurement  ac- 
curacy of  all  the  sensors.  The  WARF  data  set  is  not  limited  to  the 
results  presented  in  this  paper.  Other  analyses  of  the  radar  data  that 
were  not  obtained  in  the  vicinity  of  the  buoy  are  also  available.  These 


94 


% 

experiments  indicated  thnt  during  a hurricane,  HF  skywave  radar  can  pro- 
vide operational  surface  data  that  are  as  accurate  as  the  more  recog- 
nized in  situ  measurements.  The  supportive  surface  data  supplied  by 
y the  WARF  radar  would  prove  particularly  useful  for  tracking  during 

early  formative  stages  of  hurricanes  when  multiple  centers  may  be  ob- 
served or  when  cirrus  shielding  may  obscure  visual  location  by  satel- 
lite cloud  photography.  The  high  resolution,  large  coverage  area, 
real-time  steering, and  continuous  monitoring  capabilities  are  unique  to 
skywave  radar.  The  hurricane  data  obtained  from  skywave  radar  comple- 
ments data  obtained  from  satellites,  aircraft,  and  buoys. 


Table  3 

BABE  WIND  SPEED  ESTIMATES  MADE  ON  4 SEPTEMBER  1977 


Sensor 

Observation 

(GMT) 

Bearing 

From 

Hurricane 

Center 

(°N) 

Radial 

Distance 

From 

Hurricane 

Center 

(km) 

Wind 

Speed 

(m/s) 

WARF  (10  m) 

2253 

340 

60 

18.2 

WARF  (19.5  m) 

2253 

340 

60 

18.3 

Aircraft  (305  m) 

1800 

25 

60 

30.8 

Aircraft  (19.5  m) 

1800 

25 

60 

18.5 

Oil  Platform  (19.5  m) 

2100 

320 

100 

15.4 

Oil  Platform  (19.5  m) 

2100 

0 

35 

11.3 

WCM3  (9‘ 1 

- 

- 

60 

16.4 

WCSb  (9.1  m) 

- 

- 

60 

16.8 

Computed  maximum  wind  speed  from  Eq.  (3-35)  in  Reference  30. 

**  Computed  maximum  sustained  wind  speed  from  Eq.  (3-34)  in  Reference  30. 
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Joseph  W.  Maresca,  Jr. 

SRI  International 

ABSTRACT 

The  capability  of  HF  skywave  radar  to  measure  surface  winds,  waves, 
and  currents  at  distances  up  to  3000  km  for  tropical  and  extratropical 
storm  conditions  is  summarized.  Significant  wave  height  and  wave  spec- 
tral estimates  made  using  the  Wide  Aperture  Research  Facility  skywave 
radar  were  compared  to  similar  measurements  obtained  by  National  Data 
Buoy  Office  buoys  EB20  (41°N,  138°W)  and  EB71  (26°N,  93.5°W).  Agree- 
ment to  within  10%  was  found  for  wave  heights  under  varying  conditions, 
ranging  from  less  than  1 m to  hurricane  wave  conditions  greater  than  5 m. 


I INTRODUCTION 

High-frequency  (HF)  skywave  radar  estimates  of  the  surface  wind 
speed  and  direction,  surface  current,  significant  wave  height,  and  wave 
spectrum  can  be  made  over  several  million  square  kilometers  of  ocean  by 
remotely  measuring  the  Doppler  spectrum  of  the  sea-echo  signal.  Measure- 
ments made  at  the  SRI-operated  Wide  Aperture  Research  Facility  (WARF)  HF 
skywave  radar  demonstrate  that  the  skywave  radar  can  track  large  extra- 
tropical  storms,  and  small  intense  tropical  storms  and  hurricanes;  the 
radar  can  also  provide  detailed  surface  maps  of  the  winds,  waves,  and 
currents  within  a storm.  The  large  coverage  area  of  the  WARF  skywave 
radar  is  possible  because  the  HF  radio  waves  are  transmitted  to,  and 
returned  from,  ocean  areas  by  means  of  one  or  more  ionospheric  "reflec- 
tions." For  one  ionospheric  reflection,  ocean  areas  up  to  3000  km  away 
from  the  radar  can  be  monitored.  The  purpose  of  this  paper  is  to  sum- 
marize the  capabilities  of  HF  skywave  radar  that  utilizes  ionospheric 
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reflection,  with  special  emphasis  on  the  measurement  of  significant  wave 
height  and  the  one-dimensional  wave-frequency  spectrum. 

The  skywave  radar  estimates  of  the  surface  wave  parameters  are  ob- 
tained from  the  HF  sea-echo  Doppler  spectrum.  An  example  of  the  sea- 
echo  spectrum  recorded  at  WARF  for  51.2  s of  coherent  integration  is 
shown  in  Figure  1.  The  spectrum  consists  of  two  strong  first-order 
echoes  produced  by  the  resonant  interaction  of  ocean  waves  of  wave  num- 
ber k and  radio  waves  of  wave  number  kQ.  For  near  grazing  angles, 

k = 2k  . 

o 

In  addition,  a second-order  continuum  that  is  sensitive  to  changes  in 
the  directional  ocean  wave  spectrum  surrounds  the  first-order  echoes. 
Barrick  (1972a, b)  has  derived  theoretical  expressions  that  accurately 
describe  the  HF  scattering  process.  These  expressions  have  been  mathe- 
matically inverted  to  recover  parts  of  the  directional  wave  spectrum 
directly  from  the  Doppler  spectrum. 

Barrick  and  Lipa  (1978)  summarize  the  theory  of  HF  scattering  from 
ocean  waves;  the  methods  for  computing  the  wave  height,  the  wave  spec- 
trum, and  the  surface  current;  and  the  results  of  several  experiments 
conducted  to  validate  these  methods  through  the  use  of  HF  groundwave 
radar.  In  this  paper,  experiments  conducted  at  the  WARF  skywave  radar 
to  determine  the  accuracy  of  these  methods,  including  a simple  method 
of  mapping  wind  direction,  are  described.  The  WARF  radar  measurements 
have  been  compared  to  in  situ  point  measurements  made  from  National 
Data  Buoy  Office  (NDBO)  data  buoys  or  research  vessels.  The  relative 
agreement  between  the  radar  and  in  situ  measurements  have  been  found  to 
be  within  the  measurement  accuracy  of  the  in  situ  instrument. 


II  WARF  SKYWAVE  RADAR 

WARF  is  a bistatic  HF  skywave  radar  located  in  central  California. 
The  radar  is  operated  in  the  HF  band  between  6 and  30  MHz.  A 20-kW 
swept-frequency  continuous-wave  (SFCW)  signal  is  transmitted  from  the 
transmitter  site  at  Lost  Hills,  California;  after  round-trip  ionospheric 
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FIGURE  1 RANGE-DOPPLER-PROCESSED  SEA-ECHO  DOPPLER  SPECTRUM.  The  mean 

Doppler  spectrum  is  an  average  of  Doppler  spectra  recorded  at  different  range  lines 
separated  by  3 km.  The  first-order  echoes  produced  by  a resonant  interaction 
between  the  radio  waves  and  the  ocean  waves  is  sensitive  to  changes  in  the  wind- 
direction  field.  The  second-order  sideband  structure  surrounding  the  stronger  Bragg 
line  is  sensitive  to  changes  in  the  directional  ocean  wave  spectrum. 


propagation,  the  signal  is  received  at  the  receiving  site  185  km  to  the 
north  at  Los  Banos,  California.  The  receiving  array  is  2.5-kra  long  and 
forms  a beam  of  0.5°  at  15  MHz. 

The  WARF  radar  coverage  area  is  shown  in  Figure  2.  The  radar  can 
be  electronically  steered  in  0.25°  increments  anywhere  within  the  cover- 
age area.  The  minimum  range  or  skip  distance  of  the  radar  is  typically 
800  km,  and  the  maximum  range  for  one  ionospheric  reflection  is  generally 
3000  km  or  less.  The  nominal  position  accuracy  is  about  20  km,  and 
relative  accuracies  between  consecutive  measurements  are  at  least  an 
order  of  magnitude  better.  The  range  to  the  ocean  path  illuminated  is 
a function  of  time  and  radar  frequency.  In  general,  ranges  out  to  2200 
km  can  be  monitored  95%  of  the  time;  coverage  beyond  2200  km  drops  to 
50%. 

The  size  of  the  ocean  patch  monitored  by  the  WARF  radar  is  a func- 
tion of  the  sampling  parameters  and  can  be  specified  by  the  radar  opera- 
tor. The  characteristic  length-scale  of  the  ocean  wave  conditions  is 
considered  in  selecting  the  size  of  the  ocean  patch  to  be  illuminated 
by  the  radar.  For  example,  the  ocean  patch  monitored  for  a hurricane 
would  be  smaller  than  the  ocean  patch  monitored  for  a Pacific  storm. 

The  minimum  resolution  cell  routinely  achievable  is  3 km  in  range  by 
15  km  in  cross-range.  The  cross-range  resolution  is  a function  of  the 
radar  range  and  beamwidth.  For  a maximum  one-hop  range  of  3000  km,  a 
0.5°  beam  would  result  in  a cross-range  distance  of  about  26  km.  Gen- 
erally, 21  resolution  cells  or  independent  Doppler  spectrum  measurements 
spaced  at  3-km  increments  make  up  the  illuminated  ocean  patch.  If  the 
Doppler  spectra  at  each  range  interval  is  averaged,  the  size  of  ocean 
patch  will  be  63  x 25  km.  For  Pacific  Ocean  measurements,  the  63  x 25  km 
scattering  patch  is  used  for  wave  measurements;  for  hurricanes,  a smaller 
ocean  patch,  15  x 25  km,  is  used. 

Ill  SURFACE  WIND  DIRECTION  MAPS 

The  capability  and  technique  for  mapping  the  surface  wind  direction 
field  for  large  weather  systems  over  the  ocean  has  previously  been 
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demonstrated  (Long  and  Trizna,  1973;  Barnum  et  al.,  1977).  This  mapping 
technique  has  also  been  applied  directly  to  tropical  storm  and  hurricane 
wind  fields  (Maresca  and  Barnum,  1978;  Maresca  and  Carlson,  1977). 

Wind  direction  is  estimated  from  the  power  ratio  of  the  first- 
order  echo  returns  (Stewart  and  Barnum,  1975;  Barnum  et  al.,  1977). 
Agreement  between  the  WARF  radar  and  the  anemometer  measurements  of 
wind  direction  over  the  Pacific  Ocean  is  +16°  (Stewart  and  Barnum,  1975). 
For  hurricane  winds,  agreement  between  the  WARF  radar  and  NDBO  buoy  wind 
direction  measurements  is  better  than  10°  (Maresca  and  Carlson,  1977b). 
Figure  3 shows  an  example  of  the  WARF-inferred  surface  wind  map  produced 
as  Hurricane  Eloise  passed  between  NDBO  buoys  EB04  and  EB10  in  the  Gulf 
of  Mexico.  The  WARF-inferred  wind  direction  maps  have  been  used  to  lo- 
cate the  center  of  hurricanes.  The  earliest  measurements  were  made  on 
Eloise,  and  they  agreed  to  within  35  km  of  the  smooth  track  produced  by 
the  National  Hurricane  Center  (NHC) . Recently,  a skywave  radar  track 
was  compiled  for  Hurricane  Anita  from  17  wind  direction  maps  recorded 
over  a 4-day  period.  The  mean  difference  between  the  WARF  position 
estimates  and  the  NHC  smooth  track  was  19  km  (Maresca  and  Carlson, 

1978b). 

IV  SURFACE  CURRENT  CAPABILITY 

The  radial  component  of  the  surface  current  can  be  inferred  from 
the  measured  phase  velocity,  or  Doppler,  of  the  ocean  waves  producing 
the  first-order  echoes  (Barrick  et  al.,  1974;  Stewart  and  Joy,  1975; 
Barrick  et  al.,  1977).  Unlike  the  skywave  radar  wind  direction  measure- 
ment, the  effects  of  the  ionospheric  motion  must  be  known  to  make  a 
skywave  measurement  of  current.  The  entire  Doppler  spectrum  can  be 
shifted  by  ionospheric  motion.  It  has  been  shown  by  Maresca  et  al. 

(1976)  and  Maresca  and  Carlson  (1977b)  that  land  and  offshore  oil  plat- 
form echoes  received  during  coastal  scans  along  the  Gulf  of  Mexico  are 
sufficient  to  remove  the  effects  of  the  ionosphere  from  the  data.  The 
analysis  of  the  ocean  surface  would  then  be  identical  to  the  groundwave 
radar  analysis.  The  requirement  of  a reference  would  generally  limit 
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the  measurement  to  coastal  regions.  The  expected  accuracy  of  the  surface 
current  estimates  is  about  10  cm/s. 


V SIGNIFICANT  WAVE  HEIGHT  AND  WAVE  SPECTRUM  ESTIMATES 

The  measurement  of  the  significant  wave  height  and  wave  spectrum  is 
made  from  the  weak  second-order  region  of  the  sea-echo  Doppler  spectrum. 
The  accuracy  of  the  wave  estimates  is  dependent  not  only  on  the  method 
used  to  compute  the  wave  parameters,  the  number  of  radar  samples 
averaged,  and  the  natural  variability  of  the  waves,  but  also  on  the  con- 
tamination effects  of  the  ionosphere.  The  quality  of  the  skywave  radar 
data  depends  on  the  ionospheric  propagation  path.  Ionospheric  smearing 
effects  and  multiple-path  propagation  can  destroy  the  second-order 
spectral  contributions  (Maresca  and  Barnum,  1977;  Maresca  and  Carlson, 

1977) .  In  general,  the  potential  contamination  caused  by  the  ionosphere 
would  be  the  largest  source  of  error  in  the  measurement.  Recent  work  by 
SRI  and  NOAA  (Georges  and  Maresca,  1978;  and  Maresca  and  Georges,  1978) 
has  resulted  in  improved  methods  for  collecting  uncontaminated  data  and 
for  recovering  the  wave  spectra,  even  when  the  data  is  contaminated  by 
the  ionosphere. 

Each  of  the  HF  inversion  formulas  summarized  from  the  measured 
Doppler  spectrum  by  Barrick  and  Lipa  (1978)  to  compute  the  ocean  wave 
spectrum  must  be  tested  for  skywave  data  (Barrick,  1977a, b;  Lipa,  1977, 

1978) .  Presented  below  are  results  of  three  wave-height  verification 
experiments  conducted  at  the  WARF  radar.  In  one  of  these  experiments, 
the  one-dimensional  wave  frequency  spectrum  is  also  computed  from  the 
radar  data.  Significant  wave  heights  during  the  experiments  ranged 
from  less  than  1 m to  more  than  5 m.  The  radar  estimates  were  compared 
to  NDBO  data  buoy  measurements  in  the  Gulf  of  Mexico  (EB71)  and  Pacific 
Ocean  (EB20).  Agreement  between  the  NDBO  buoy  measurements  and  the  WARF 
radar  measurements  for  significant  wave  heights  were  within  the  accuracy 
of  the  NDBO  buoy  measurements. 


VI  WAVE  SPECTRAL  ESTIMATES  FROM  THE  HF  DOPPLER  SPECTRUM 


Barrick  (I977a,b)  derived  the  following  approximate  closed-form 
expression  to  compute  the  one-dimensional  wave  frequency  spectrum,  S(w), 
and  rms  wave  height,  h*,  from  the  Doppler  spectrum: 
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where  u - (uiB|v-l|)  is  the  radian  ocean  wave  frequency;  Wg  is  the  radian 
Bragg  frequency;  wD  is  the  radian  Doppler  frequency;  v - “>D/u)B;  a^WgV) 
and  (^((DgV)  are  the  first-  and  second-order  power  contribution  to  the 
Doppler  spectrum  expressed  as  radar  cross  section  per  mean  surface  area 
per  radian  per  second  of  bandwidth;  kQ  is  the  radian  radio  wave  number; 
and  W(v)  is  a weighting  function  derived  by  Barrick  (1977a) 

The  weighted  second-order  power  contribution  used  in  Eq.  (1)  is 
obtained  by  averaging  the  return  from  both  sides  of  the  stronger  first- 
order  echo.  The  rms  wave  height  is  obtained  by  dividing  the  total 
weighted  second-order  power  surrounding  both  sides  of  the  stronger 
first-order  echo  by  the  total  first-order  power.  If  one  assumes  that 
the  process  is  Gaussian,  the  significant  wave  height,  H is  4h*. 

The  accuracy  of  these  expressions  can  be  evaluated  by  inverting 
theoretical  Doppler  spectra  prodqced  from  known  input  directional  wave 
spectra.  The  results  indicate  that  Eqs.  (1)  and  (2)  tend  to  over- 
estimate wave  height  by  a constant  amount.  Barrick  (1977a)  has  shown 
that  when  the  ratio  of  the  actual  wave  height  to  the  radar-measured  wave 
height  h/hfc,  is  plotted  as  a function  of  the  parameter  k^^,  h/hfc  is  a 
constant  for  kQh  > 0.20.  Maresca  and  Carlson  (1977b)  also  showed  the 
dependence  of  h/h*  on  radar  frequency  and  wave  directionality. 
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Barrick  (1977a)  has  tested  Eq.  (2)  by  using  surface  wave  HF  radar 
data  and  in  situ  buoy  measurements,  and  the  agreement  was  found  to  be 
within  22%.  The  rms  error  is  primarily  dependent  on  the  accuracy  of 
the  method  used  to  compute  wave  height  and  the  number  of  independent 
Doppler  spectra  averaged  before  computing  the  wave  height.  Each  Doppler 
spectrum  of  the  sea-echo  signal  in  Barrick' s test  was  an  average  of  9 
samples.  If  the  number  of  samples  averaged  had  been  increased,  then 
better  estimates  of  the  sea  echo  would  have  resulted  and  the  22%  error 
might  have  been  reduced. 


VII  WARF  MEASUREMENTS  OF  SIGNIFICANT  WAVE  HEIGHT  AND  THE  WAVE  SPECTRUM 

The  results  of  three  WARF  wave  height  verification  experiments  con- 
ducted on  2 June  1977,  23  August  1977,  and  1 September  1977  are  reported 
here  to  demonstrate  the  accuracy  of  the  skywave  radar  measurement  tech- 
nique over  a range  of  seas  and  ionospheric  conditions.  The  radar  mea- 
surements were  recorded  in  the  vicinity  of  NDBO  buoys  EB71  (26°N,  93.5°W) 
in  the  Gulf  of  Mexico  and  EB20  (41°N,  138°W)  in  the  Pacific  Ocean.  EB71 
is  located  approximately  2900  km  from  WARF;  EB20  is  located  about  1700  km 
from  WARF.  The  Doppler  spectra  were  produced  from  102.4  s of  coherent 
integration.  Eq.  (1)  was  used  to  compute  the  ocean  wave  spectrum,  and 
Eq.  (2)  was  used  to  compute  the  rms  wave  height.  The  significant  wave 
height  was  computed  by  multiplying  the  rms  wave  height  by  4.0.  At 
least  40  Independent  samples  of  the  Doppler  spectra  were  averaged.  The 
errors  in  the  spectral  estimates  are  proportional  to  1//N  , where  N is 
the  number  of  samples  (Barrick  and  Snider,  1977).  The  sampling  error 
based  on  40  or  more  averages  for  all  three  experiments  is  16%  or  less. 

The  nominal  accuracy  of  NDBO  moored  buoys  measurement  of  significant  wave 
height  reported  by  the  Data  Quality  Division  of  NDBO  is  ±0.3  m.  The 
WARF  measurements  of  significant  wave  height  made  for  all  three  experi- 
ments agree  to  within  ±0.3  m of  the  buoy  measurements.  The  results  of 
some  NDBO  calibration  experiments  are  described  in  Steele  and  Johnson 
(1978). 
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A.  June  2,  1977  Experiment 

The  placement  of  NDBO  data  buoy  EB20  on  station  in  the  Pacific 
Ocean  in  1977  provided  the  first  opportunity  to  compare  the  WARF- 
estimates  of  significant  wave  height,  Hg,  with  an  independent  in  situ 
measurement.  The  Doppler  spectra  were  taken  on  2 June  1977  to  the  west 
of  an  atmospheric  front.  The  front  passed  through  EB20,  travelling  to 
the  east.  During  a period  of  less  than  6 hours,  Hg  increased  from  1.4  m 
to  3.0  m and  remained  constant  at  approximately  3 m for  more  than  6 
hours.  The  significant  wave  height  measured  over  a 12-hour  period  is 
given  in  Table  1 . 

The  WARF  radar  measurement  reported  here  was  centered  125  km  west 

of  EB20  at  1819Z.  Since  the  wave  height  and  wave  spectral  measurements 

at  EB20  remained  constant  (±0.1  m)  from  1800  to  2400Z,  we  expected  that 

the  WARF  measurement  of  the  significant  wave  height  should  be  about  3.0  m. 

We  compared  the  WARF  radar  estimates  of  significant  height  computed  from 

Eq.  (2)  and  the  one-dimensional  wave  frequency  spectrum  computed  from 

Eq.  (1)  to  those  measurements  made  at  EB20.  The  WARF  estimate  of  H is 

s 

3.0  m;  the  WARF  measured  wave  spectrum  is  shown  in  Figure  4.  We  com- 
pared these  estimates  to  those  made  at  EB20  at  2100Z  because  the  WARF 
measurement  was  located  west  of  the  buoy  where  the  seas  were  probably 
slightly  larger.  However,  the  difference  in  the  EB20  wave  spectra 
measurements  at  1800Z  and  2100Z  are  negligible,  and  comparison  of  the 
EB20  and  WARF  measurements  any  time  between  1800Z  and  2100Z  produces 
essentially  identical  results.  The  agreement  between  the  significant 
wave  heights  measured  at  EB20  and  WARF  is  within  the  ±0.3  m accuracy  of 
the  NDBO  moored  buoys.  Further  details  of  this  experiment  are  available 
in  Maresca  and  Carlson  (1977a)  and  Maresca  and  Georges  (1978b). 

B.  August  23,  1977  Experiment 

The  placement  of  the  NDBO  data  buoy  EB71  on  station  in  1977  provided 
the  first  opportunity  to  compare  WARF  estimates  of  significant  wave 
height  in  the  Gulf  of  Mexico  with  an  in  situ  measurement.  WARF  measure- 
ments recorded  at  2140Z  on  23  August  1977  were  centered  on  EB71.  Wave 
conditions  at  EB71  remained  fairly  constant  before,  during,  and  after 
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Table  1 


SIGNIFICANT  WAVE  HEIGHT 

MEASURED  AT  EB20  JUNE  2,  1977 

* 

Time 

Significant  Wave 

(GMT) 

Height  (m) 

1200 

1.4 

1500 

1.8 

- 

1800 

3.0 

2100 

3.1 

2400 

2.9 

Table  2 

SIGNIFICANT  WAVE  HEIGHT 

' 

MEASURED  AT  F.B71 

AUGUST  23,  1977 

Time 

Significant  Wave 

(GMT) 

Height  (m) 

1200 

0.9 

1500 

0.9 

1800 

1.0 

2100 

0.9 

2400 

0.8 

] 


'1 

I ; 

i 
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the  experiment,  as  shown  in  Table  2.  Wave  conditions  were  generally 
1.0  m or  less.  The  WARF  estimate  of  Hg  is  1.0  m,  which  compares  reason- 
ably well  with  the  0.9  m measurement  made  at  2100Z  at  EB71.  Theoretical 
simulations  suggest  th  .t  the  accuracy  of  the  method  used  to  calculate 
h*  decreases  significantly  for  Hg  <1.0  m,  and  this  experiment  probably 
represents  a lower  limit  on  the  measurement  capability. 

C.  September  1,  1977  Experiment 

Hurricane  Anita  was  tracked  by  the  WARF  radar  from  29  August 

through  2 September  1977.  WARF  measurements  of  H were  made  at  2324Z 

s 

on  31  August  1977  in  the  right  rear  quadrant  of  Anita  as  the  storm 
passed  near  EB71.  The  WARF  measurement  was  centered  about  80  km  from 
the  radar -derived  storm  center.  Comparison  of  the  WARF  estimates  of  H 

s 

with  EB71  estimates  of  H was  difficult  because  small  differences  in 

s 

location  of  the  measurement  with  respect  to  the  center  of  the  hurricane, 

where  conditions  change  rapidly  over  short  distances,  can  result  in 

different  values  of  Hg.  The  areal  distribution  of  the  significant  wave- 

height  field  was  generated  from  EB71  measurements  of  Hg  as  the  storm 

passed  the  buoy.  It  was  assumed  that  the  spatial  distribution  of  the 

hurricane  wave  conditions  remained  constant  for  ±12  hours.  The  WARF 

estimate  of  significant  wave  height  was  then  compared  to  an  interpolated 

value  of  Hg  computed  from  the  EB71  significant  wave  height  field  (Figure 

5).  The  WARF  estimate  of  Hg  was  5.4  m,  and  the  interpolated  value  of 

H from  the  EB71  estimates  was  5.5  m.  This  is  within  the  0.3  m error 
s 

associated  with  the  NDBO  buoy  measurement.  Estimates  of  the  WARF 
significant  wave  heights  computed  in  all  four  quadrants  of  the  hurricane 
will  be  described  in  more  detail  in  Maresca  and  Carlson  (1978a). 


VIII  SUMMARY  AND  CONCLUSIONS 

The  main  purpose  of  this  paper  was  to  summarize  the  capability  of 
measuring  significant  wave  height  and  the  wave  spectrum  by  HF  skywave 
radar.  Agreement  between  the  radar  estimates  of  significant  wave 
height  and  NDBO  data  buoy  measurements  was  within  10%  for  a wide  range 
of  wave  conditions.  Similar  accuracy  was  found  for  the  wave  spectrum 
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FIGURE  5 COMPARISON  OF  THE  SIGNIFICANT  WAVE  HEIGHT  MEASURED  BY  THE 
WARF  SKYWAVE  RADAR  AND  NDBO  EB71.  The  EB71  significant  wave 
height  estimate*  were  compiled  from  24  hours  of  data  and  were  plotted  with 
respect  to  the  hurricane  center. 
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estimate.  The  demonstrated  capability  by  skywave  radar  for  continuously 
tracking  and  monitoring  the  surface  winds  and  waves  throughout  all  re- 
gions of  a hurricane  is  unique.  One  skywave  radar  similar  in  design  to 
WARF  could  routinely  monitor  the  incident  wave  conditions  along  entire 
continental  shelf  regions  off  the  east  or  west  coast  of  the  United 
States. 

ACKNOWLEDGMENTS 

This  work  was  funded  in  part  by  the  Air  Force  Office  of  Scientific 
Research  (AFOSR)  and  in  part  by  the  NOAA  Wave  Propagation  Laboratory. 

The  wave  spectral  data  was  obtained  from  the  National  Data  Buoy  Office. 


References 


Barnum,  J.  R. . , J.  W.  Maresca,  Jr.,  and  S.  M.  Serebreny,  1977.  High- 

resolution  mapping  of  oceanic  wind  fields  with  skywave  radar,  IEEE 
Trans,  on  Antennas  and  Propagation,  AP-25,  128-132. 

Barrick,  D.  E.,  1972a.  First-order  theory  and  analysis  of  MF/HF/VHF 

scatter  from  the  sea,  IEEE  Trans,  on  Antennas  and  Propagation,  AP- 

20,  2-10. 

Barrick,  D.  E.,  1972b.  Remote  sensing  of  sea  state  by  radar,  in:  Remote 
Sensing  of  the  Troposphere,  (ed.  V.  E.  Derr),  Chapter  12,  U.S. 
Government  Printing  Office. 

Barrick,  D.  E.,  J.  M.  Headrick,  R.  W.  Bogle,  and  D.  D.  Crombie,  1974. 

Sea  backscatter  at  HF:  interpretation  and  utilization  of  the  echo. 
Proceedings  of  the  IEEE,  62,  673-680. 

Barrick,  D.  E.  and  J.  B.  Snider,  1977.  The  statistics  of  HF  sea-echo 
Doppler  spectra,  IEEE  Trans,  on  Antennas  and  Propagation,  AP-25, 
19-28. 

Barrick,  D.  E. , M.  W.  Evans,  and  B.  L.  Weber,  1977.  Ocean  surface  cur- 
rents mapped  by  radar.  Science,  198,  138-144. 

Barrick,  D.  E.,  1977a.  Extraction  of  wave  parameters  from  measured  HF 
sea-echo  Doppler  spectra,  Radio  Science,  12,  415-424. 

Barrick,  D.  E.,  1977b.  The  ocean  waveheight  nondirectional  spectrum 

from  inversion  of  the  HF  sea-echo  Doppler  spectrum.  Remote  Sensing 
of  Environment,  6,  201-227. 

Barrick,  D.  E.  and  B.  J.  Lipa,  1978.  Ocean  surface  features  observed  by 
HF  coastal  ground-wave  radars:  a progress  review.  Ocean  Wave  Climate, 
(eds.  M.  D.  Earle  and  A.  Malahoff),  Plenum  Publishing  Company, 

New  York. 

Georges,  T.  M.  and  J.  W.  Maresca,  Jr.,  1978.  The  effect  of  radar  beam- 
width  on  the  quality  of  sea-echo  Doppler  spectra  measured  with  HF 
skywave  radar,  in  preparation. 

Lipa,  B.  J.,  1977.  Derivation  of  directional  ocean-wave  spectra  by 
integral  inversion  of  second-order  radar  echoes.  Radio  Science, 

12,  425-434. 

Lipa,  B.  J.,  1978.  Inversion  of  second  order  radar  echoes  from  the  sea, 

J.  Geophys.  Research,  83,  959-962. 

Long,  A.  E.  and  D.  B.  Trlzna,  1973.  Mapping  of  North  Atlantic  winds  by 
HF  radar  sea  backscatter  interpretation,  IEEE  Trans,  on  Antennas 
and  Propagation,  AP-21,  680-685. 


117 


Maresca,  J.  W.  Jr.,  J.  R.  Barnum,  and  K.  L.  Ford,  1976.  HF  skywave  radar 
measurements  of  coastal  and  open  ocean  surface  currents  (Abstract), 
1976  Annual  Meeting  of  USNC/URSI  University  of  Massachusetts, 
Amherst,  Massachusetts,  October  11-15,  1976. 

Maresca,  J.  W.  Jr.  and  J.  R.  Barnum,  1977.  Measurement  of  oceanic  wind 
speed  from  HF  sea  scatter  by  skywave  radar,  IEEE  Trans,  on  Antennas 
and  Propagation,  AP-25,  132-136. 

Maresca,  J.  W.  Jr.  and  C.  T.  Carlson,  1977a.  Remote  measurement  of  the 
ocean  wave  frequency  spectrum  by  HF  skywave  radar  (Abstract),  1977 
Fall  Meeting  American  Geophysical  Union,  San  Francisco,  California, 
December  5-9,  1977. 

Maresca,  J.  W.  Jr.  and  C.  T.  Carlson,  1977b.  Tracking  and  monitoring  of 
hurricanes  by  HF  skywave  radar  over  the  Gulf  of  Mexico,  Technical 
Report  1,  SRI  International,  Menlo  Park,  California. 

Maresca,  J.  W.  Jr.  and  J.  R.  Barnum,  1978.  Remote  measurements  of  the 
position  and  surface  circulation  of  hurricane  Eloise  by  skywave 
radar.  Monthly  Weather  Review,  in  press. 

Maresca,  J.  W.  Jr.  and  C.  T.  Carlson,  1978a.  HF  skywave  radar  measure- 
ments of  wave  height  for  hurricane  Anita,  in  preparation. 

Maresca,  J.  W.  Jr.  and  C.  T.  Carlson,  1978b.  HF  skywave  radar  track  of 
hurricane  Anita,  in  preparation. 

Maresca,  J.  W.  Jr.  and  T.  M.  Georges,  1978a.  Estimating  significant 
wave  height  from  the  HF  Doppler  spectrum  of  the  sea  echo  received 
from  two  or  more  ionospheric  paths,  in  preparation. 

Maresca,  J.  W.  Jr.  and  T.  M.  Georges,  1978b.  HF  skywave  radar  measure- 
ment of  the  ocean  wave  spectrum,  in  preparation. 

Maresca,  J.  W.  Jr.  and  T.  M.  Georges,  1978c.  Ionospheric  diagnostics 
and  sea  scatter  propagation  management,  in  preparation. 

Steele,  K.  and  A.  Johnson,  Jr.,  1978.  NDBO  wave  measurements.  Ocean 
Wave  Climate,  (eds.  M.  D.  Earle  and  A.  Malahoff),  Plenum  Pub- 
lishing Corporation,  New  York. 

Stewart,  R.  H.  and  J.  R.  Barnum,  1975.  Radio  meapurements  of  oceanic 
winds  at  long  ranges:  an  evaluation.  Radio  :ience,  10,  853-857. 

Stewart,  R.  H.  and  J.  W.  Joy,  1974.  HF  radio  measurements  of  surface 
currents,  Deep-Sea  Research,  21,  1039-1049. 


*w 


118 


